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SECTION  1 
INTRODUCTION 

The  present  Miss  Distance  Indicator  (MIDI)  Radar  was  developed  in  1973  by 
the  U.  S.  Army  Air  Defense  Board  and  the  Uestinghouse  Defense  and  Space 
Center.  MIDI  performs  an  essential  role  in  the  development,  test,  and 
evaluation  of  air  defense  weapons.  The  MIDI  radar  tracks  a  target  and 
measures  the  vector  miss  distance  (both  distance  and  direction)  to  the  point 
of  closest  approach  of  projectiles  fired  at  the  target,  as  illustrated  in 
Figure  1-1.  The  MIDI  is  capable  of  performing  other  range  instrumentation 
radar  functions  such  as  tracking  bullets  or  missiles  to  provide  trajectory 
data.  Two  Interchangeable  antennas  (3  foot  diameter  and  6  foot  diameter) 
permit  either  short  range,  wide  f ield-of-view  operation  or  long  range,  narrow 
f ield-of-view  operation.  This  feature  permits  scoring  of  a  wide  range  of 
projectile  types.  The  mobility  of  the  MIDI  allows  it  to  be  used  wherever 
needed  within  the  Ft.  Bliss  or  White  Sands  Missile  Range  military 
reservations. 

The  recent  air  defense  artillery  systems  have  been  developed  to  counter 
the  various  low  altitude  threats  such  as  nap-of-the-earth  penetration  aircraft 
and  attack  helicopters.  The  U.S.  Army  Air  Defense  Artillery  Board  needs  a 
Miss  Distance  Indicator  Radar  that  provides  tracking  and  scoring  of  air 
defense  weapons  against  simulated  low  altitude  threats.  Ideally  this  upgraded 
radar  should  retain  all  the  desirable  characteristics  of  the  current  system, 
e.g. ,  accuracy  and  mobility,  but  be  capable  of  scoring  more  closely  to  the 
horizon. 

Two  types  of  interference  —  clutter  and  multipath  —  degrade  the 
tracking  accuracy  against  low  altitude  targets.  Clutter  is  the  backscatter  of 
the  radar  energy  from  the  ground.  Multipath  interference  is  a  forward 
scattering  phenomenon  in  which  the  radar  echo  from  the  target  is  reflected 
from  the  ground,  creating  a  false  "image"  in  addition  to  the  direct  return  to 
the  radar.  Clutter  will  be  strongest  and  multipath  reflection  interference 
will  be  weakest  from  rough  surfaces.  Conversely,  for  smooth  terrain,  clutter 
will  be  weak  and  multipath  interference  will  dominate.  The  terrain  roughness 
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from  a  radar  standpoint  depends  not  only  upon  the  physical  features  but  also 
upon  the  radar  wavelength  and  grazing  angle.  All  surfaces  appear  smoother 
when  viewed  at  small  grazing  angles. 

The  U.S.  Army  Air  Defense  Board  contracted,  via  the  Training  and  Doctrine 
Command  Contracting  Activity  at  Ft.  Gustis,  with  Georgia  Tech  for  a 

Methodology  Improvement  Program  to  examine  techniques  and  methodologies  that 
will  improve  low  angle  elevation  tracking  and  scoring  capabilities  of  a  second 
generation  MIDI  radar.  The  major  output  of  this  program,  as  stated  in  the 
contract,  is  a  report  stating: 

1.  Conclusions  and  recommendations  concerning  the  best  method  or 

combination  of  methods  to  employ  for  accurate  low-elevation  tracking 
and  scoring  and 

2.  Minimum  achievable  elevation  angle  specified  in  terms  of 

technological  risks  as  a  function  of  elevation  angle. 

The  general  approach  was  as  follows: 

1.  Become  familiar  with  current  and  future  MIDI  requirements.  This  was 

achieved  through  study  of  several  documents  provided  by  the  U.S.  Army 
Air  Defense  Board,  visits  by  Army  personnel  to  Georgia  Tech, 
observing  MIDI  missions  at  the  North  McGregor  Test  Range,  and 

numerous  discussions  with  Ft.  Bliss  and  White  Sands  Missile  Range 
personnel. 

2.  Eliminate  the  less  desirable  techniques  based  upon  MIDI  operational 
requirements. 

3.  Calculate  the  performance  for  the  most  promising  methods  by  varying 
several  radar  parameters  (e.g.,  antenna  beamwidth,  pattern  shape, 
number  of  Doppler  filters),  geometric  parameters  (e.g.,  radar  height 
and  target  height)  and  environmental  parameters  (e.g.,  backscatter 
coefficient  and  reflection  coefficient). 

4.  Measure  the  environmental  parameters  at  the  N:th  McGregor  Test  range 
using  MIDI  to  the  extent  possible  within  program  constraints  such  as 
range  time  available  and  weather. 


5.  Validate  the  hardware  feasibility  of  achieving  the  required 
transmitter  power,  digital  signal  processing  throughput  rates,  and 
storage. 

The  conclusions  and  recommendations  were  presented  at  a  program  review  in 
a  letter  report  during  August  1984.  This  final  report  reiterates  these 
results  and  provides  detailed  descriptions  of  the  supporting  analyses. 


SECTION  2 

MAJOR  RADAR  PERFORMANCE  OBJECTIVES 


The  desired  MIDI  radar  operating  modes  are  summarized  In  Table  2-1.  The 
primary  use  of  the  system  Is  tracking  the  unmanned  targets  and  measuring 
vector  miss  distance  during  live  firing  weapon  tests.  Air  defense  weapon 
scoring  Is  desired  within  the  range/altitude  envelopes  shown  in  Figures  2-1 
and  2-2.  No  augmentation,  passive  or  active.  Is  permitted  on  the  bullets  or 
missiles.  The  radar  must  be  capable  of  tracking  an  Individual  bullet  or 
missile  to  provide  trajectory  data.  Another  use  is  tracking  manned  aircraft 
for  mission  control  or  range  safety.  The  manned  aircraft  and  the  recoverable 
unmanned  targets  can  be  equipped  with  a  radar  beacon  to  improve  detection  and 
tracking  in  clutter.  The  expendable  towed  targets  will  be  passively  augmented 
to  a  radar  cross  section  of  one  square  meter  or  greater. 

Gun  1 

Gun  1  fires  20  millimeter  projectiles  at  rates  of  3600  per  minute.  The 
projectile  velocity  at  the  target  vicinity  may  be  as  high  as  1000  meters  per 
second.  Firings  of  intermixed  types  of  ammunition,  e.g. ,  tracer  rounds,  with 
slightly  different  ballistics  can  cause  the  bullets  to  bunch  or  pass.  The 
MIDI  should  score  all  ammunition  loads  properly.  The  required  scoring 
accuracy  for  Gun  1  is  one  milliradian.  This  translates  to  a  miss  distance 
accuracy  of  1  meter  at  a  target  range  of  1000  meters. 

Gun  2 

Gun  2  fires  40  millimeter  projectiles  at  rates  of  less  than  1000  per 
minute.  The  other  requirements  are  the  same  as  those  of  Gun  1. 


MisSiles 


The  performance  envelope  boundaries  for  missile  scoring  were  established 
for  the  long  range,  vehicular-mounted  launcher  class  of  air  defense 
missiles.  A  missile  diameter  of  15  inches  was  used  in  estimating  the  radar 
cross  section.  The  required  angular  scoring  accuracy  for  missiles  scoring,  is 
0.1  milliradian;  this  translates  to  a  miss  distance  accuracy  of  1  meter  at  a 
target  range  of  10  kilometers. 


Figure  2-2  Gun  Scoring 


TABLE  2-1.  DESIRED  MIDI  OPERATING  MODES 


OBJECT 

BEACON 

TRACKING 

MAXIMUM 

ENVELOPE 

MINIMUM 

TYPICAL 
SCORING  RANGES 

TRACKED 

MODE 

RANGE 

ALTITUDE 

BULLET 

MISSILE 

BY  RADAR 

(km) 

(m) 

(km) 

(km) 

UNMANNED  TARGETS 

TOW 

30 

150 

1-4 

2-20 

HELICOPTER 

30 

15 

1-4 

2-20 

FIXED  WING 

YES 

30 

15 

1-4 

2-20 

MANNED  AIRCRAFT  YES  30 


15 


BULLETS 


6 


MISSILES 


30 


SECTION  3 
CLUTTER  ANALYSIS 


The  present  noncoherent  MIDI  radar  system  Is  designed  to  operate  at 
elevation  angles  of  3.5  beamwidths  (BW)  and  above.  The  system  employs  no 
moving  target  indication  (MTI)  processing  for  clutter  interference 
reduction.  Instead,  the  required  clutter  rejection  is  achieved  by  attenuation 
of  the  received  clutter  power  through  the  reduced  gain  of  the  antenna 
elevation  sidelobes.  At  the  minimum  3.5  BW  elevation  angle,  the  antenna 
provides  a  two-way  attenuation  of  greater  than  50  dB  which  is  adequate  to 
remove  clutter  interference  as  a  serious  source  of  tracking  error. 

The  second  generation  MIDI  must  operate  at  elevation  angles  considerably 
lower  than  3.5  BW.  Elevation  angles  as  low  as  0.2°  are  desirable.  At  these 
low  angles,  the  received  clutter  power  will  be  entering  the  radar  through  the 
main  beam  of  the  antenna,  thus  removing  the  antenna  pattern  as  a  mechanism  for 
clutter  rejection.  Consequently,  the  required  clutter  rejection  must  be 
achieved  by  signal  processing  in  the  receiver.  Because  of  the  large  velocity 
differences  between  the  targets  of  interest  and  the  clutter  background,  MTI 
Doppler  processing  provides  an  excellent  alternative  for  clutter  reduction. 
The  following  section  provides  a  detailed  analysis  of  the  nature  of  the  MIDI 
clutter  background,  the  types  of  MTI  processing  that  could  be  employed  to 
reduce  clutter  interference,  and  the  effects  of  the  residual  clutter  on  the 
second  generation  MIDI  scoring  performance. 

3.1  MIDI  CLUTTER  REFLECTIVITY  MEASUREMENTS 

During  a  visit  to  the  present  MIDI  radar  location  at  the  North  McGregor 
Missile  Range,  New  Mexico,  Georgia  Tech  personnel  utilized  the  MIDI  system  to 
collect  samples  of  ground  clutter  video  data  at  four  ranges  and  three 
elevation  angles.  Because  of  instrumentation  vans  located  in  the  vicinity  of 
the  MIDI  system,  only  one  azimuthal  direction  provided  an  unobstructed  view  of 
the  desert  surface.  The  resulting  set  of  twelve  data  runs  is,  then  lore,  by 
no  means  a  definitive  description  of  ground  clutter  at  the  MIDI  site,  but  it 
proved  useful  in  determining  reasonable  bounds  for  clutter  reflectivity 
coefficients  in  a  representative  MIDI  environment. 


3.1.1  CLUTTER  MEASUREMENT  AND  ANALYSIS  PROCEDURE 


The  clutter  measurement  procedure  was  designed  to  produce  sum  channel 
video  clutter  data  as  well  as  sufficient  calibration  information  to  allow 
reduction  of  the  raw  clutter  data  to  equivalent  clutter  reflectivity  (m^/m^) 
information.  The  calibration  procedure  consisted  of  two  stages.  First,  the 
amplitude  response  function  of  the  receiver  was  determined  before  and  after 
data  collection  using  10  dB  step  attenuators  built  into  the  MIDI  system.  The 
response  of  the  reciver  from  noise  level  to  70  dB  above  noise  was  determined 
by  boresighting  the  antenna  on  a  calibration  beacon  and  collecting  video 
samples  as  receiver  attenuation  was  varied  from  0  to  70  dB.  This  documented 
the  receiver  response  over  the  signal  range  of  interest.  Second,  a  reference 
sample  of  video  data  was  collected  while  the  antenna  was  boresighted  on  a 
trihedral  corner  reflector  located  on  a  pole  6  km  north  of  the  MIDI.  This 
provided  a  sample  of  output  power  for  a  known  radar  cross  section  (RCS)  target 
at  a  known  distance. 

Once  the  calibration  was  completed,  the  radar  antenna  was  pointed  in  the 
obstruction-free  direction  and  clutter  data  collection  was  begun.  At  each 
elevation  position,  the  range  gate  was  set  to  various  ranges  having  a  large 
clutter  return  and  several  seconds  of  video  data  were  collected.  Therefore, 
the  samples  collected  represented  a  higher  clutter  power  than  would  have 
resulted  from  sampling  every  range  bin.  The  collected  data  thus  represent  a 
worst  case  condition.  This  consideration  was  included  in  the  analysis  of  the 
required  clutter  rejection. 

The  reduction  of  the  clutter  data  to  useful  reflectivity  values  involved 
three  steps.  Calibration  of  the  received  video  clutter  data  with  the  receiver 
response  function  resulted  in  clutter  power  values  in  dB  above  noise. 
Normalizing  the  clutter  power  data  using  the  range  and  the  reference  corner 
reflector  video  resulted  in  clutter  RCS  (a)  in  dBsm.  Finally,  the  clutter 

RCS  was  normalized  with  respect  to  illuminated  area  to  produce  clutter 

O  2  2 

reflectivity  (o  )  values  in  dB  (m  /m  ). 


3.1.2  CLUTTER  MEASUREMENT  RESULTS 


The  results  of  the  MIDI  clutter  measurements  are  presented  in  Table  3- 

l.  Four  different  ranges  were  sampled,  varying  from  ~  1,500  m  to  -  25,000 

m.  Three  elevation  settings  corresponding  to  horizontal,  one  beamwidth  above 
horizontal,  and  one  beamwidth  below  horizontal  were  used  at  each  range.  As 
can  be  seen  from  the  table,  the  reflectivity  values  do  not  follow  any  regular 
pattern.  Raising  and  lowering  the  antenna  at  a  given  range  resulted  in  as 
much  as  14  dB  of  variation  in  the  return  signal.  This  was  a  result  of  the 
uneven  nature  of  the  terrain  surface,  particularly  at  the  longer  ranges,  where 
the  terrain  became  hilly.  In  fact,  the  high  reflectivity  values  at  25  km  are 
discounted  because  they  represent  backscatter  data  from  the  foothills  of  a 
local  mountain  range.  The  rest  of  the  data  exhibit  a  reflectivity  that  varies 
from  -15  dB  to  -38  dB.  The  following  clutter  interference  analysis  uses 
values  for  o°  of  -20  dB  and  -30  dB  which  adequately  represents  the  expected 
average  values. 

3.2  MIDI  SCORING  PROJECTILE  RCS 

The  second  generation  MIDI  must  perform  scoring  missions  for  a  large 
variety  of  weapons,  ranging  from  the  20  mm  Vulcan  gun  up  to  the  Patriot  class 
missile  at  ranges  from  0.5  km  to  30  km.  Four  basic  weapon  types  were 
considered:  a  20  mm  gun,  a  40  mm  gun,  and  two  air  defense  missiles.  A 
detailed  tracking  error  analysis  requires  estimation  of  the  radar  cross 
section  of  these  projectiles. 

The  MIDI  Final  Technical  Report  produced  by  Westinghouse  reports  the  RCS 
of  a  20  mm  projectile  to  be  -30  dBsm. ^  The  source  of  this  value  was  not 
referenced,  but  since  it  was  used  in  designing  the  original  MIDI,  the  same 
value  was  used  in  our  analysis.  A  means  of  predicting  the  X-band  RCS  of  the 
40  mm  projectile  was  developed  by  comparing  the  nominal  -30  dBsm  RCS  of  the  20 
mm  projectile  with  theoretical  predictions  as  follows.  The  RCS  Handbook 
edited  by  George  Ruck  was  examined  for  theoretical  and  experimental  RCS  values 
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"Mobile  Miss  Distance  Indicator  (MIDI)  Final  Technical  Report," 
Westinghouse  Electric  Corporation,  Data  Item  A003,  Contract  DAAD07-72-C- 
0155,  May  1973. 


Range 

Elevation 

RCS 

RCS 

Reflectivity 

Cm) 

(mils) 

Cm2) 

(dBsm) 

(m2/m2) 

(dB) 

for  objects  of  the  same  general  shape  as  bullet  projectiles.^  A  bullet  is 
most  closely  modeled  as  a  flattened  ogive;  however,  the  RCS  Handbook  presents 
no  data  for  such  an  object.  A  short  cylinder  was  considered  the  next  best 
approximation,  but  no  RCS  data  existed  for  a  cylinder  with  both  dimensions  in 
the  resonant  region,  which  represents  20  mm  and  40  mm  bullets  at  X-band. 
Consequently,  the  initial  approach  was  to  estimate  the  RCS  from  a  circular 
flat  plate  at  normal  incidence.  This  was  used  to  represent  the  RCS  of  the 
bullet  projectile  when  viewed  from  the  rear. 

Figure  7-24  from  the  RCS  Handbook  presents  radar  cross  sections  for  a 
disk  at  normal  Incidence  as  a  function  of  wave  number  kQ  and  radius  a.  This 
graph  is  reproduced  in  Figure  3-1.  Table  3-2  demonstrates  the  results  of  this 
RCS  calculation  for  the  two  projectiles  of  interest: 


TABLE  3-2.  RCS  OF  A  DISK  AT  NORMAL  INCIDENCE 


Projectile 

a  (m) 

V 

o  (M2) 

a  (dBsm) 

20  mm  bullet 

0.01 

2.0 

0.003 

-25 

40  mm  bullet 

0.02 

4.0 

0.02 

-17 

The  theoretical  estimate  of  -25  dBsm  for  the  20  mm  projectile  is  5  dB  higher 
than  the  RCS  used  in  the  MIDI  design.  Presumably,  this  5  dB  accounts  for  the 
reduction  in  RCS  at  aspect  angles  off  of  normal.  Applying  the  same  correction 
factor  to  the  theoretical  40  mm '  projectile  RCS  of  -17  dBsm  results  in  an 
average  cross  section  of  -22  dBsm.  The  values  of  -30  dBsm  and  -22  dBsm  are 
used  in  the  tracking  error  analyses  that  follow. 

The  radar  cross  sections  of  potential  missiles  to  be  scored  by  the  MIDI 
was  less  easily  derived.  No  data  for  RCS  measurements  of  advanced  air  defense 
missiles  was  found  by  Georgia  Tech,  and  any  experimental  RCS  value  would 
almost  certainly  be  classified.  Because  of  the  complex  shape  of  missiles. 
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Ruck,  G.  T. ,  Editor  Radar  Cross  Section  Handbook,  Plenum  Press,  1970 


Source:  Ruck  (1970) 


Figure  3-1.  Exact  and  physical-optics  backscatter  cross 
sections  for  a  perfectly  conducting  disc  at 
normal  incidence. 

Target  Velocity  (m/s) 


15  30 


T/rget 

Doppler  Frequency 
(Hz) 


Figure  3-2.  Improvement  factor  of  a  Butterworth  delay 
line  canceller. 


including  Che  concave  nozzle  of  the  propulsion  unit,  a  simple  but  accurate 
theoretical  prediction  using  physical  optics  was  not  possible.  The  source 
used  for  missile  RCS  estimates  was  the  MIDI  Final  Technical  Report ^  This 

report  stated  that  the  RCS  of  missiles  could  range  from  -17  dBsm  to  -5  dBsm  at 
X-band.  The  larger  RCS  value  of  -5  dBsm  was  used  to  represent  the  long  range 
air  defense  missiles. 

3.3  CLUTTER  REJECTION  REQUIREMENTS 

The  representative  clutter  reflectivity  values  of  -20  dB  and  -30  dB  were 
used  to  predict  the  amount  of  clutter  power  received  by  the  MIDI  radar  under 
various  configurations.  The  expected  projectile  RCS  and  the  required  signal- 
to-clutter  ratio  (SCR)  were  used  to  estimate  the  necessary  MTI  improvement 
factor. 

The  slgnal-to-interference  ratio  (SIR)  required  for  adequate  performance 
of  the  MIDI  tracking  radar  was  derived  using  procedures  described  in  the 
Handbook  of  Radar  Measurements  by  Barton  and  Ward^.  In  Chapter  8,  Barton 
outlines  a  method  for  calculating  the  required  SNR  using  the  desired  accuracy 
and  the  radar  system  parameters.  The  same  procedure  can  be  applied  to  the 
(SCR)  ratio  as  follows. 

The  required  energy  ratio,  R,  is  given  by 

0,  2 

R  -  (~J->  ,  3.1 

,  Qm 

where 

0^  ■  the  3  dB  beamwidth, 

km  -  the  normalized  monopulse  slope,  and 

o„  ■  the  desired  rras  tracking  error. 

O 

From  the  energy  ratio  the  SIR  is  computed  from  the  equation 


3.  Barton,  D.  K.  and  Ward,  H.  R. ,  Handbook  of  Radar  Measurements,  Prentiss- 
Hall,  1969. 


where 


n 


the  number  of  pulses  on  the  target  and 
the  filter  mismatch  loss. 


The  3  dB  beamwldth  of  the  present  MIDI  large  antenna  Is  22  mrad,  and  the 
monopulse  error  slope  Is  approximately  1.8.  Instrumentation  and  mechanical 
Inaccuracies  limit  the  minimum  rms  tracking  error  to  approximately  0. 1  mrad. 
For  these  values.  Equation  3.1  produces  a  required  energy  ratio  of  15,000. 
The  number  of  pulses  on  a  projectile  such  as  a  bullet,  is  limited  to 
approximately  20  by  the  projectile  velocity.  If  the  receiver  bandwidth  B  is 
matched  to  the  pulse  length  t  such  that  Bt  *  1  ,  then  the  filter  mismatch 
loss  will  be  1  dB.  For  these  values,  Equation  3.2  predicts  a  required  SIR  of 
26  dB. 

The  MTI  improvement  factor  required  to  produce  a  26  dB  SCR  is  given  by 

x  ,  SCRolA  3.3 

°T 

where 

A  *  R  sin  03  •  -  area  illuminated, 

o°  ■  clutter  reflectivity, 

aT  *  target  RCS,  and 

SCR  -  required  signal-to-clutter  ratio. 

Table  3-3  demonstrates  the  required  improvement  factor  for  a  variety  of 
scenarios.  The  MTI  improvement  values  range  from  a  low  of  36  dB  for  missile 
scoring  at  30  km  in  low  clutter  to  a  high  of  60  dB  for  20  mm  bullet  scoring  at 
2  km  in  high  clutter.  An  improvement  factor  of  60  dB  is  difficult  to  achieve, 
particularly  with  as  few  pulses  as  are  available  to  the  MIDI  system.  For  this 
reason,  the  MTI  processing  analysis  that  follows  concentrates  on  achieving 
improvement  factors  of  40  and  50  dB. 
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Table  3-3  REQUIRED  MTI  IMPROVEMENT  FACTOR  FOR  A  26  dB  SCR 


Projectile 

RCS 

(dBsm) 

Range 

(km) 

BW 

(deg) 

0 

a 

(dB) 

I 

(dB) 

20  mm  bullet 

-30 

2 

1.25 

-20 

60 

-30 

50 

40  mm  bullet 

-22 

4 

1.25 

-20 

54 

-30 

44 

Missile 

-5 

30 

1.25 

-20 

46 

-30 

36 

3.4  MTI  PROCESSING  ANALYSIS 

The  high  radial  velocity  of  missiles  and  bullets  to  be  scored  by  the  MIDI 
and  the  stationary  arrangement  of  the  radar  make  an  ideal  situation  for 
implementation  of  Doppler  processing.  The  very  speed  that  makes  MTI 
processing  attractive,  however,  induces  technical  problems  that  impact  the 
selection  and  design  of  candidate  MTI  systems.  This  section  discusses  the 
nature  of  the  Doppler  detection  problem,  the  candidate  MTI  techniques 
available,  and  the  types  of  improvement  factors  that  can  be  achieved. 

3.4.1  PRF  SELECTION 

The  projectiles  that  the  second  generation  MIDI  must  score  may  be 
travelling  at  speeds  up  to  1000  m/s.  This  velocity  corresponds  to  a  required 
unambiguous  Doppler  frequency  f^  of  64,500  Hz  at  X-band.  Some  form  of  PRF 
staggering  must  be  employed  to  extend  the  unambiguous  Doppler  frequency  f^  to 
the  required  range  of  64,500  Hz.  For  a  two  PRF  stagger,  the  appropriate  pulse 


The  maximum  scoring  projectile  velocity  is  1000  m/s.  The  sampling  time 
tQ  is  determined  by  the  time  the  projectile  is  within  the  5  m  range  gate. 
Since  the  range  gates  are  tracking  the  drone  target,  which  may  have  velocities 
on  the  order  of  200  m/s,  the  net  velocity  difference  between  projectile  and 
range  gate  could  be  as  great  at  1200  m/s.  For  a  5  m  range  gate,  this 
correspond  to  a  sampling  time  of  4.17  ms.  The  number  of  pulses  on  the 
projectile  would  then  be 

N  -  PRF  .  tQ  3.4 

*  PRF  •  4.17  ms. 

For  a  two  PRF  stagger,  the  maximum  unambiguous  frequency  occurs  when 
multiples  of  each  PRF  overlap.  This  occurs  at  a  frequency  fj  given  by 

prfl  •  prf2 

f d  *  PRF  "  -  PRF~  *  3,5 

By  establishing  a  known  relationship  between  PRFj  and  PRF2,  Equations  3.4  and 
3.5  can  be  used  to  determine  the  appropriate  PRF  values.  If  a  bank  of  N 
narrow  band  filters  were  employed  to  provide  velocity  information,  the  two 
PRF's  would  be  chosen  such  that  the  first  ambiguous  frequency  of  PRFj  would 

fall  in  the  middle  of  the  N  -  1  filter  of  PRF2.  This  statement  corresponds  to 

the  following  equation: 

PRFj  -  prf2.  3.6 

A  quadratic  equation  for  PRF2  is  obtained  by  combining  equations  3.4, 
3.5,  and  3.6 

t  PRF^  -  PRF.  -  f  -  0 
o  l  l  a 


The  solution  of  Equation  3.7  is: 


If  f0  4.17  ms  and  f<j  ■  64,500  Hz,  Equation  3.8  gives  PRF2  -  4055  Hz.  Then  by 
Equation  3.5,  PRF^  =  3815.  For  a  1000  m/s  projectile,  these  PRF's  will 
provide  on  the  order  of  16  pulses  on  targets.  This  is  a  worst  case 
estimate.  The  average  bullet  velocity  in  the  scoring  region  is  approximately 
680  m/s.  At  this  velocity,  32  pulses  could  be  accumulated.  This  result  is 
considered  in  the  following  MTX  analysis. 

3.4.2  MTI  PROCESSING  TECHNIQUES 

Two  MTI  processing  techniques  considered  in  this  analysis  can  provide  the 
necessary  1000  m/s  coverage  in  the  Doppler  domain:  a  delay  line  canceller  and 
a  coherent  FFT.  The  characteristics  of  these  filters  are  described  in  the 
following  sections,  and  their  relative  merits  as  applied  to  the  second 
generation  MIDI  radar  are  considered. 

3.4.2. 1  Delay  Line  Canceller 

Delay  line  cancellers  (DLCs)  are  a  common  technique  for  achieving  a  broad 
passband  for  detecting  moving  targets  while  filtering  out  a  large  proportion 
of  the  stationary  clutter  return.  Since  the  DLC  can  be  implemented  digitally, 
a  large  number  of  range  gates  can  be  accommodated,  and  good  MTI  improvement 
can  be  achieved  with  three  or  four  delay  stages.  Increasing  in  the  number  of 
stages  increases  the  length  of  time  necessary  for  the  canceller  to  achieve  a 
steady-state  condition  and  places  more  stringent  requirements  on  transmitter 
and  oscillator  stability.  But  increasing  the  number  of  stages  results  in  a 
sharper  slope  of  the  MTI  response  function  and  thus  increases  the  clutter 
cancellation  produced  by  the  filter.  The  response  function  H(f)  of  a 
Butterworth  canceller  is  closely  approximated  by  the  equation 

| H(f ) | 2  =  f  2m  3.9 

1  +  (tt) 

where  fc  =  the  cutoff  frequency  and 

m  -  the  number  of  poles  or  stages  in  the  filter. 
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The  improvement  factors  for  3-stage  and  a  4-stage  Butterworth  DLC  as  a 
function  of  target  velocity  is  presented  in  Figure  3-2.  Both  cancellers 
achieve  their  maximum  improvement  at  a  target  velocity  of  15  m/s.  The 
improvement  of  the  3-stage  canceller  is  limited  to  40  dB,  but  the  4-stage 
filter  has  a  theoretical  limit  of  50  dB.  The  4-stage  canceller  is 
considerably  more  effective  in  detecting  small  bullet  targets  in  high 
background  clutter. 

The  DLC  filters  discussed  so  far  are  for  a  f ixed  PRF  and  have  a  passband 
limited  to  the  PRF.  For  the  second  generation  MIDI  to  fulfill  its  scoring 
mission  on  high  speed  projectiles,  pulse-to-pulse  staggering  of  the  PRF  must 
be  employed  to  extend  the  blind  velocity  Vg  to  the  1000  m/s  requirement.  At 
an  average  PRF  of  4000  Hz  the  first  unstaggered  blind  speed  V  occurs  at 

v  rl  +  r2  +  r3  +  ...  +  rn 


where  Rj ,  R2  . . . ,  Rjj  are  the  ratios  of  the  pulse  intervals  expressed  as 


integers1 


A  typical  velocity  response  curve  for  a  staggered  DLC  from 


Skolnik's  Radar  Handbook  is  presented  in  Figure  3-3. The  overlap  of  blind 
velocities  of  the  individual  PRF’s  causes  modulation  of  the  response  function 
in  the  passband  region.  The  severity  of  this  modulation  is  reduced  as  the 
number  of  PRF’s  employed  increases.  For  the  second  generation  MIDI,  the  ratio 
of  true  blind  speed  to  the  unstaggered  blind  speed  is  V/Vg  ■ 
"'6'2^m7s'8  *  adequate  Doppler  coverage  could  be  achieved  with  a  pulse 

internal  ratio  of  19:16:17:13. 

Staggering  the  pulse  interval  has  drawbacks  as  well  as  advantages.  In  a 
nonstaggered  multistage  DLC,  a  slowly  varying  waveform,  such  as  a  specular 
clutter  return,  will  be  almost  perfectly  cancelled.  But  in  a  staggered  DLC, 
the  time  interval  between  the  two  samples  will  vary,  resulting  in  a  voltage 
residue  out  of  the  canceller.  Skolnik  has  derived  the  MTI  improvement  factor 


where 


n  *  number  of  pulses  on  target  and 
6  »  ratio  of  pulse  intervals. ^ 

For  the  staggering  choice  previously  described  the  maximum  improvement  factor 
would  be  41  dB.  If  this  improvement  were  deemed  insufficient,  a  5  interval 
PRF  stagger  might  be  required. 
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Source:  Skolnik  (1970) 
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Figure  3-3.  Velocity- response  curve:  dual  canceller,  no 
feedback,  11:16:13:17  pulse-interval  ratio. 


The  previous  analysis  indicates  that  a  4-stage  delay  line  canceller  (DLC) 
with  4  or  5  period  PRF  stagger  could  achieve  the  improvement  factors  required 
in  the  MIDI  II  scenario.  Delay  line  cancellers  have  a  disadvantage  which 
would  seriously  limit  their  usefulness  in  the  MIDI  scenario.  As  previously 
mentioned,  DLC  filters  require  a  period  of  several  pulses  to  achieve  their 
optimum  steady  state  condition.  This  effect  is  demonstrated  in  Figure  3-4 
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which  presents  the  response  of  a  DLC  to  a  transient  signal.  This  response  is 
representative  of  a  bullet  passing  through  a  scoring  range  gate  at  600  m/s. 
The  canceller  takes  4  or  5  pulses  to  reach  optimum  operation,  and  another  4  or 
5  pulses  to  return  to  zero  once  the  bullet  is  gone.  This  degradation  of 
performance  would  be  acceptable  if  it  were  the  only  source  of  transient 
errors.  Since  the  scoring  range  gates  of  the  MIDI  system  would  track  the 
drone  target  at  speeds  up  to  200  m/s,  large  specular  clutter  sources  could 
suddenly  appear  within  a  range  gate  even  if  the  objects  themselves  were 
stationary.  Figure  3-5  depicts  the  canceller  response  to  such  a  transient 
clutter  signal.  While  the  DLC  effectively  removes  the  clutter  during  steady 
state  operation,  the  sudden  appearance  and  disapperance  of  the  clutter  signal 
results  in  strong  spikes  in  the  output  voltage.  Under  certain  conditions, 
these  transient  responses  to  clutter  could  greatly  increase  the  number  of 
clutter  false  alarms.  This  would  seriously  affect  the  MIDI  scoring 
performance.  Consequently,  the  delay  line  canceller  is  not  recommended  as  a 
solution  to  clutter  interference  in  the  second  generation  MIDI  radar  system. 


3.4. 2. 2  Coherent  Fast  Fourier  Transform 

Recent  advances  in  the  speed  of  digital  microprocessors  and  other 

hardware  have  made  fast  Fourier  transforms  (FFT)  increasingly  useful  in  radar 
applications  where  high  data  rates  are  an  inherent  requirement.  Like  the  DLC, 
the  FFT  can  be  implemented  digitally,  allowing  high  speed  performance  at  a 
moderate  cost.  The  FFT,  however,  possess  a  number  of  characteristics  that 
make  it  superior  to  a  delay  line  canceller. 

One  PRF  interval  of  the  coherent  video  spectrum  of  a  moving  target  in 
noise  and  clutter  is  depicted  in  Figure  3-6.  The  bandwidth  of  the  target 

spectrum  (Bs)  is  primarily  determined  by  the  length  of  the  sampling  period 
(Bg  a  l/tQ),  whereas  the  bandwidth  of  the  clutter  spectrum  is  determined  by 
both  the  sampling  time  and  internal  clutter  motion.  For  the  MIDI  radar,  the 
average  sampling  time  is  so  short  (~  5  msec)  that  it  will  almost  certainly  be 
the  limiting  factor  on  clutter  spectrum  bandwidth  as  well.  The  fast  Fourier 
transform  essentially  acts  as  a  bank  of  Np  contiguous,  narrowband  Doppler 

filters  ranging  from  -PRF/2  to  +PRF/2.  Optimally,  the  bandwidth  of  these 


filters  is  chosen  to  equal  the  target  spectrum  bandwidth  Bs;  this  choice 
maximizes  the  signal-to-noise  ratio  (SNR)  out  of  each  filter.  In  fact,  one  of 
the  advantages  of  the  FFT  is  that  the  SNR  out  of  the  filter  is  the  product  of 
the  SNR  at  the  input  to  the  filter  and  the  number  of  filters.  The  delay  line 
canceller,  on  the  other  hand,  passes  nearly  all  of  the  noise  through  its  wide 
passband.  An  additional  advantage  of  the  filter  bank  is  that  the  target's 
velocity  can  be  determined  from  the  filter  in  which  it  was  detected.  This 
could  be  useful  in  matching  projectile  trajectories  to  precomputed  ballistic 
tables. 

The  FFT  has  a  final  advantage  over  the  DLC.  Because  it  is  a  batch 
process,  rather  than  a  continuous  one,  the  FFT  does  not  display  the  adverse 
response  to  transient  clutter  and  target  signals  produced  by  the  DLC.  The 
result  of  a  sudden  appearance  of  a  strong  clutter  signal  during  a  batch  will 
be  a  widening  of  the  clutter  spectrum  associated  with  the  reduced  sampling 
time.  Such  a  widening  could  reduce  target  detectability  In  the  two  filters 
contiguous  with  the  clutter  filter,  but  it  will  not  have  the  disastrous  effect 
produced  by  the  DLC  circuit. 

Coherent  FFT  Improvement  Factor 

The  improvement  factor  I  of  the  FFT  can  be  computed  as  follows: 


where  and  C^  are  the  target  and  clutter  powers  into  the  FFT  processor  and 
SQ  and  CQ  are  the  respective  powers  out  of  the  processor.  The  FFT  filter  bin 
size  was  assumed  to  be  matched  to  the  target  bandwidth.  Signal  losses  due  to 
filter  mismatch  as  well  as  to  other  sources  in  the  processor  are  considered 
later. 

The  frequency  response  of  the  FFT  for  a  single  frequency  bin  is  pictured 
in  Figure  3-7.  The  uniform  sidelobes  are  a  result  of  a  Dolph-Chebyshev 
weighting  applied  to  the  input  time  samples.  Some  form  of  weighting  is 
necessary  in  any  practical  FFT  application  to  reduce  the  frequency  sidelobes; 
Dolph-Chebyshev  weighting  is  commonly  employed  because  it  produces  uniform 
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Figure  3-6.  One  PRF  interval  of  the  coherent  video  spectrum  of 
a  moving  target  in  noise  and  clutter. 
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Figure  3-7.  Filter  response  of  target  frequency  bin  with  Dolph-Chebyshev 
weighting  applied. 


height  sidelobes.  The  FFT  response  was  evaluated  for  two  frequency  sldelobe 
levels:  -35  dB  and  -50  dB.  The  amount  of  clutter  power  out  of  the  FFT 
processor  is  the  sum  of  the  clutter  power  within  the  frequency  bin  of  the 
target  plus  the  residual  clutter  that  is  passed  through  the  sidelobes  of  the 
filter  response. 

The  improvement  factor  produced  by  a  coherent  FFT  is  demonstrated  in 
Figure  3-8  for  a  16-point  and  a  32-point  FFT.  In  both  cases,  the  maximum 
improvement  achieved  is  limited  by  the  frequency  sidelobe  level  (FSL).  Since 
adequate  MIDI  II  tracking  performance  requires  40  dB  to  50  dB  of  improvement, 
the  frequency  sidelobes  should  be  suppressed  as  much  as  possible.  As  can  be 
seen  in  the  figure,  the  narrower  frequency  bins  of  the  32-point  FFT  result  in 
a  more  rapid  rise  to  the  maximum  improvement  factor.  In  adddition,  they  allow 
a  more  accurate  estimation  of  the  target  velocity,  which  could  prove  useful  in 
determining  vector  miss  distance. 

Coherent  FFT  Losses 

There  are  a  number  of  sources  of  signal  loss  in  the  implementation  of  a 
fast  Fourier  Transform.  The  process  of  weighting  the  input  signal  train  to 
reduce  the  frequency  sldelobe  level  (FSL)  results  in  a  loss  in  signal  power. 
The  lower  the  desired  sldelobe  level,  the  higher  the  loss  associated  with 
it.  For  the  Dolph-Chebyshev  weighting  described  previously,  the  loss  is 
0.8  dB  for  35  dB  sidelobes  or  1.4  dB  for  50  dB  sidelobes. 

In  addition  to  the  weighting  loss,  there  is  a  filter  mismatch  loss  that 
results  from  target  energy  extending  beyond  the  bandwidth  of  the  FFT  bin  in 
which  it  resides.  The  amount  of  this  loss  is  a  function  of  the  FFT  bin  size 
and  the  width  of  the  target  spectrum.  The  FFT  bin  sizes  considered  in  this 
analysis  were  250  Hz  for  the  16-point  FFT  and  125  Hz  for  the  32-point  FFT. 
For  the  high  speed  targets  of  the  MIDI  system,  the  target  spectrum  bandwidth 
will  be  determined  by  the  sampling  time  tQ;  i.e. 


This  bandwidth  could  be  as  high  as  240  Hz  for  1000  m/s  missiles  shot  at  an 
oncoming  drone  target.  A  more  representative  bandwidth  would  be  120  Hz,  which 
corresponds  to  the  average  velocity  of  bullets  (600  m/s)  in  the  vicinity  of 
the  drone  aircraft.  For  a  120  Hz  target  spectrum  bandwidth,  the  16-point  FFT 
has  a  filter  mismatch  loss  of  1.5  dB;  i.e.,  30%  of  the  power  extends  beyond 
the  bandwidth  of  the  FFT  bin.  For  the  the  same  target  signal,  the  32-point 
FFT  has  a  filter  mismatch  loss  of  4  dB.  In  this  case,  more  than  half  of  the 
energy  extends  beyond  the  filter.  These  loss  figures  assume  that  the 
detection  threshold  is  applied  to  each  bin  individually.  Some  of  the  lost 
energy  could  be  recovered  by  integrating  over  several  bins  surrounding  the 
target  bin.  Barring  such  processing,  the  total  loss  associated  with  the  FFT 
signal  processing  is  summarized  in  Table  3-4. 

TABLE  3-4.  COHERENT  FFT  PROCESSING  LOSSES 


Source 

16-point  FFT 

(dB) 

32-point 

(dB) 

FFT 

Mismatch  loss 

1.5 

1.5 

4.0 

-35  dB  FSL 

-50  dB  FSL 

0.8 

1.4 

0.8 

1.4 

Total  losses 

2.3 

2.9 

4.8 

5.4 

In  spite  of  the  processing  losses  associated  with  FFT  processing,  it  is 
still  the  best  approach  for  the  second  generation  MIDI  design.  The  desired  50 
dB  improvement  factor  can  be  achieved  with  proper  weighting  to  reduce 
frequency  sidelobes.  The  1000  m/s  unambiguous  velocity  can  be  reached  with  a 


two  stage,  batch  staggered  PRF  of  4055  and  3815  Hz.  Because  the  FFT  is  a 
batch  process,  the  undesirable  transient  response  of  clutter  signals  which 
affects  the  DLC  design  is  absent.  For  all  these  reasons,  a  coherent  FFT  is 
the  best  MTI  processing  design  for  the  MIDI  II. 

3.5  CLUTTER  INTERFERENCE  TRACKING  ERROR  ANALYSIS 

The  low  angle  scoring  coverage  desired  for  the  second  generation  MIDI 
radar  places  stringent  requirements  on  the  MTI  processor.  A  graphic  depiction 
of  the  coverage  requirements  were  presented  in  Figure  2-1.  The  effect  of 
residual  clutter  interference  on  MIDI  tracking  performance  was  analyzed  for 
these  low  angle  scenarios.  The  results  were  compared  to  the  desired  coverage 
diagram  to  determine  the  effectiveness  of  MTI  processing  in  reaching  the  low 
angle  goals. 

3.5.1  CLUTTER  INTERFERENCE  MODEL 

The  effect  of  clutter  interference  on  MIDI  II  tracking  performance  was 
simulated  with  a  computer  model.  The  physical  scenario  for  the  model  is  shown 
in  Figure  3-9.  A  drone  target  is  flown  toward  the  radar  at  a  constant 
altitude  ht.  The  radar  sits  on  a  pedestal  at  a  height  of  5  m  above  the 
surrounding  terrain.  The  radar  is  assumed  to  be  boresighted  on  the  drone 
target,  and  the  scoring  target  (projectile  or  missile)  flies  past  the  drone  at 
some  offset  angle  0T«  The  clutter  power  returned  from  the  same  range  cell  as 
the  scoring  target  enters  the  antenna  at  an  angle  0^,  .  The  gains  associated 
with  these  signals  are  listed  in  Table  3-5. 


TABLE  3-5.  TARGET  AND  CLUTTER  ANTENNA  GAINS 


The  clutter  RCS  is  computed  from  the  reflectivity  o  and  the  Illuminated 
area  A: 


0£  ■  o°  A  .  3.13 

For  the  gains  previously  described,  the  received  signal-to-clutter  ratios 
(SCR)  in  the  sum  and  difference  channels  were  computed  as  follows: 


and 


(-) 

Vzin 


GZ  °T 

Gz  °C 


S  Gz  (0^) 

'c^in  ■  Oj  <0C)  <ec)  oc  • 


3.14 


3.15 


where  oT  *  scoring  target  RCS.  Finally,  the  MTI  processing  in  the  receiver 
provides  an  increase  in  the  SCR  out  of  the  signal  processor  such  that 
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The  modified  signal-to-clutter  ratios  were  used  in  the  following  computer 
program  to  compute  the  tracking  error  introduced  by  the  clutter  interference. 


3.5.2  MTI  ANGLE  PROCESSING 

Unlike  conventional  monopulse  radars,  which  process  sum  and  difference 
channels  separately,  the  MIDI  monopulse  system  converts  the  basic  A 
and  Z  channels  into  (Z/2  +  jA  /2  )  and  (Z/2  -j  A//2  )  channels,  as 
illustrated  in  Figure  4-3.  The  phase  angle  between  the  two  channels  is 
computed  from 


♦x  <9) 


2  arctan 
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where  A(0)  ”  difference  voltage  at  angle  9 


Z(9)  ■  sum  voltage  at  angle  0,  and 
L  »  loss  difference  between  A  and  Z  channels. 

The  MIDI  system  uses  a  look-up  table  to  compute  the  of f-boreslght  angle  0^ 
from  the  measured  phase  angle 

3.5.3  EFFECT  OF  CLUTTER  INTERFERENCE  ON  MIDI  ANGLE  PROCESSING 

Signals  received  from  clutter  or  multipath  will  interfere  with  an 
accurate  estimate  of  the  of f-boresight  angle  location  of  scoring  targets. 
Note  that  an  estimate  of  target  location  will  still  be  made  in  the  presence  of 
interference,  but  the  result  will  be  erroneous.  Since  the  radar  gives  no 
indication  that  interference  is  present,  it  is  important  that  the  useful 
limits  of  the  MIDI  II  system  be  determined. 

Interfering  signals  from  clutter  backscatter  will  contaminate  the 

true  Z  and  A  target  signals  as  depicted  in  Figure  3-10.  The  amplitude  of  the 

clutter  signal  will  depend  on  the  radar  cross  of  the  clutter  and  the  clutter 

attenuation  of  the  signal  processor.  The  phase  of  the  clutter  signal  relative 

to  the  target  signal,  6  ,  will  vary  randomly  from  0  to  360°.  The  contaminated 

c 

voltage  signals  will  thus  be 

I'  -  ZT  +  Zc  cos  6C  3.18 

and 

A'  *  At  +  Ac  cos  6C  ,  3.19 

target  sum  channel  voltage, 
target  difference  channel  voltage, 
clutter  sum  channel  voltage, 
clutter  difference  channel  voltage,  and 
relative  clutter-to-target  phase  angle. 

These  signals  can  be  related  to  the  pure  target  signal  and  the  clutter-to- 
slgnal  voltage  (C/S)v  ratio  by 

l'  ”  Zm( 1  +  cos  6„) 


where  ZT  - 
*T  " 
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Figure  3-9.  Physical  scenario  for  clutter  interference  model, 
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Figure  3-10.  Effect  of  clutter  interference  on  MIDI  angle  processing. 
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The  contaminated  phase  angle  measurement  Is  then 

/2  A_  1  +  (r-)][  cos  5. 
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The  MIDI  look-up  table  Is  used  to  translate  the  phase  angle  A  to  a  geometric 

E 

angle  9„.  The  measured  angle  error  Induced  by  clutter  Interference  is  then 
E 

A8  -  0E  -  eT.  3.23 

3.5.4  COMPUTER  PROGRAM  "CLTERR" 

A  computer  model  was  generated  to  compute  angle  error  as  a  function  of 
range  for  a  varitey  of  target  and  environmental  characteristics.  Table  3-6 
lists  the  input  variables  used  to  generate  the  matrix  of  performance 
predictions  presented  in  Appendix  D.  A  clutter  reflectivity  of  -20  dB  was 
used  for  all  cases.  This  represents  a  worst  case  condition  for  the  present 
MIDI  environment. 


in  detail  in  Section  3.4.6.  The  complete  set  of  output  plots  is  Included  in 
Appendix  D. 

3.5.5  CLUTTER  TRACKING  ERROR  RESULTS 

The  output  of  the  clutter  interference  tracking  error  model  consisted  of 
plots  of  rms  tracking  error  in  mils  as  a  function  of  range.  Figure  3-11 
depicts  an  example  of  this  output  for  the  case  of  an  X-band  MIDI  with  40  dB 
MTI  improvement  scoring  a  20  mm  projectile  fired  at  a  drone  target  at  50  m 
altitude.  As  can  be  seen  in  the  figure,  good  tracking  accuracy  is  achieved 
out  to  a  range  of  1700  m,  where  the  rms  error  is  less  than  0.5  mils.  Beyond 
this  range,  the  tracking  error  increases  exponentially  and  exceeds  1  mils  at  a 
range  of  2000  m.  The  sharp  rise  in  tracking  error  is  the  result  of 
illumination  of  the  clutter  area  by  the  difference  antenna  pattern.  While  the 
low  gain  portion  of  the  pattern  is  on  the  ground,  the  effect  of  clutter 
interference  changes  slowly.  As  the  illuminating  point  approaches  the  3  dB 
point  on  the  difference  pattern,  the  influence  of  clutter  increases  rapidly. 

A  significant  improvement  in  tracking  error  performance  can  be  achieved 
by  increasing  the  MTI  improvement  to  50  dB,  as  shown  in  Figure  3-12.  The 
shape  of  the  curve  is  unchanged,  but  because  of  the  increased  improvement,  the 
rms  tracking  error  at  2000  m  is  less  than  0.35  mils. 

The  tracking  performance  for  each  weapon  was  summarized  as  a  function  of 
target  height  and  range  to  present  the  results  of  the  tracking  error  analysis 
in  a  form  more  useful  to  a  radar  system  performance  evaluation.  A  total  rms 
error  of  1  mil  was  chosen  as  an  appropriate  design  specification.  For 
statistically  independent  error  sources,  this  error  allowance  would  be 
equivalent  to  four  errors  of  0.5  mil  each.  This  limit  of  0.5  mil  was  applied 
to  both  the  multipath  and  clutter  interference  analyses. 

Figure  3-13  presents  the  MIDI  II  20  mm  bullet  scoring  performance  for  a 
tracking  error  of  0.5  mil.  The  performance  limits  for  a  second  generation 
MIDI  with  40  dB  and  5(  dB  of  MTI  improvement  are  compared  to  the  coverage 
goal.  For  comparison,  the  operating  limit  of  the  present  MIDI  system  is  also 
shown.  Clearly,  the  very  low  angle  coverage  of  0.4°  is  not  achievab1  ven 
with  50  dB  of  MTI  improvement.  The  addition  of  MTI  processing,  however,  does 
provide  a  great  improvement  over  the  present  MIDI  capabilities. 
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Figure  3-13.  X-band  MIDI  20  mm  bullet  scoring:  X-band; 

40  dB  MTI  improvement;  50  m  target  height. 


TABLE  3-6.  CLUTTER  INTERFERENCE  MODEL 
TEST  MATRIX 


VARIABLE 


VALUES 


Projectile  type/RCS  (dBsm) 
Projectile  scoring  angle  (BW) 
Target  Height  (m) 

Target  Range  (km) 

MTI  Improvement  Factor  (dB) 
Clutter  Reflectivity  (dB) 


-30,  -22,  -5 
.25,  .5 
15,  50,  150 
.5  to  30 
40,  50 
-20 


The  computer  model  reproduced  the  angle  error  calculations  discussed  in 
Section  3.5.3.  The  program  computations  proceeded  as  follows: 

1.  Read  in  antenna  pattern  and  compute  MIDI  calibration  curves. 

2.  Read  in  user  selected  test  variables  (target  RCS,  height,  etc.). 

3.  Select  maximum  range. 

4.  Compute  projectile  and  clutter  antenna  gains. 

5.  Compute  signal-to-clutter  ratios  after  processing. 

6.  Compute  tracking  phase  angle  for  projectile  alone. 

7.  Compute  tracking  phase  angle  for  projectile  plus  clutter  for  200 
samples  of  clutter-to-projectile  phase  angle. 

8.  Compute  rms  tracking  error  for  uniform  clutter-to-projectile  phase 
distribution. 

9.  Output  rms  error. 

10.  Select  new  range  and  repeat  steps  4  through  9. 

11.  Plot  results. 


The  above  procedure  resulted  in  plots  of  rms  tracking  error  as  a  function 
of  range  for  a  variety  of  targets  and  conditions.  These  results  are  discussed 


A  similar  comparison  was  performed  for  the  40  mm  projectile.  Because  of 
the  increased  bullet  cross  section,  the  scoring  accuracy  achievable  at  2  km  is 
much  better  than  for  the  20  mm  projectile.  The  design  goal  for  40  mm 
projectile  scoring,  however,  was  a  15  m  altitude  target  at  4  km.  This 
corresponds  to  an  elevation  angle  of  0.2°.  As  can  be  seen  in  Figure  3-14, 
such  low  coverage  is  simply  not  possible  with  a  1.25°  beamwidth  radar  located 
behind  the  firing  line.  Some  alternative  solution  must  be  considered  for 
achieving  these  low  altitude  coverage  regions. 

The  scoring  performance  analysis  was  not  limited  to  projectile  weapons 
but  included  a  representative  missile  weapon  as  well.  The  very  low  angle 
coverage  (0.17°)  required  to  score  missiles  at  their  performance  limits 
necessitated  a  reduction  in  the  antenna  beamwidth.  An  X-band  antenna  large 
enough  to  produce  a  narrow  enough  beamwidth  would  be  physically  impractical. 
Consequently,  the  missile  scoring  performance  analysis  assumed  a  Ka~band  MIDI 
system.  A  two  meter  antenna,  such  as  in  use  in  the  present  MIDI  system,  would 
have  a  K  -band  beamwidth  of  0.35°.  This  three-fold  reduction  in  the  beamwidth 

a 

would  allow  considerably  lower  angle  tracking  coverage.  The  total  rms  error 
of  1  mil  assumed  in  the  previous  analysis,  however,  is  not  sufficiently 
restrictive  for  20  to  30  km  ranges.  At  20  km,  a  1  mil  error  would  be 
equivalent  to  20  meters.  This  is  insufficient  accuracy  for  useful  scoring 
information.  For  this  reason,  the  missile  analysis  assumed  an  rms  scoring 
error  limit  of  0.1  mil  for  clutter  interference. 

The  results  of  the  K  -band  MIDI  missile  scoring  analysis  are  presented  in 

a 

Figure  3-15  superimposed  over  the  desired  missile  coverage  diagram.  As  with 
the  projectile  weapons,  a  second  generation  MIDI  located  behind  the  firing 
line  cannot  achieve  the  low  angle  missile  coverage  desired  with  the  kind  of 
accuracy  needed  for  missile  scoring.  A  50  dB  improvement,  however,  would 
allow  accurate  scoring  for  100  meter  altitude  missile  targets  at  20  km.  This 
is  significantly  better  performance  than  could  be  achieved  with  the  present 
MIDI  system. 


N. 


SECTION  4 

MULTIPATH  INTERFERENCE  ANALYSIS 

The  primary  tasks  of  the  multipath  Interference  analysis  were: 

1.  Predict  the  MIDI  scoring  accuracy  as  a  function  of  target  altitude 
and  range,  antenna  beamwidth,  frequency,  and  terrain  surface  roughness 
(or  more  specifically,  the  ground  reflection  coefficient)  and  compare 
candidate  radar  mechanizations. 

2.  Determine  the  values  of  ground  reflection  coefficient  which 
represent  the  terrain  at  Ft.  Bliss  and  White  Sands  Missile  Range. 
Analysis  was  started  using  values  from  0.25  to  0.7.  Subsequently,  a 
computer  program  that  models  an  area  of  the  White  Sands  Missile  Range 
predicted  values  from  0.25  to  0.67.  Finally,  a  measured  value  of  0.35 
was  obtained  during  a  limited  test  program  at  North  McGregor  Range. 

3.  Determine  the  correlation  between  errors  in  target  tracking  and 
errors  in  bullet  miss  distance  measurement. 

This  section  of  the  report  describes  the  equations  used,  the  computer 
programs  that  implement  these  equations,  and  the  results. 

4. 1.  CALCULATION  OF  ANGLE  MEASUREMENT  ACCURACY 

The  MIDI  radar  determines  miss  distance  by  tracking  the  gunnery  target 
and  performing  of f-boresight  angle  measurements  of  the  projectiles  at  two 
ranges  in  the  vicinity  of  the  target. 

The  angle  measurement  errors  due  to  multipath  interference  were 
calculated  for  two  monopulse  angle  measurement  methods,  referred  to  as 
conventional  monopulse  and  MIDI  monopulse.  The  errors  were  calculated  for  two 
representations  for  interference  —  a  three  dimensional  model  providing 
specular  and  diffuse  reflections  and  a  simple  two  dimensional  specular 
reflection  model.  The  angle  measurement  models  and  the  two  computer  programs 
which  implement  the  models  are  described  herein. 
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4.1.1  ANGLE  MEASUREMENT  EQUATIONS 

4. 1.1.1  Conventional  Monopulse  The  elevation  sum  (E)  and  elevation 
difference  (A)  patterns  are  inputs  to  this  model.  The  of f-boreslght  angle 
measurement  calibration  procedure  consists  of  constructing  a  look-up  table 
of  A/E  voltage  ratios  versus  angle.  The  A/E  ratios  for  angles  below  antenna 
boresight  (negative  angles)  have  a  negative  value.  This  represents  an  antenna 
which  has  been  designed  such  that  the  response  to  a  single  signal  In  the 

A  channel  is  either  in  phase  or  180  degrees  out  of  phase  with  the 

E  response.  If  a  second  signal  such  as  multipath  interference  is  present, 
its  E  and  A  responses  are  in  phase  (or  180  degrees  out  of  phase)  with  each 
other  but  may  have  any  phase  with  respect  to  the  first  signal.  Therefore,  the 
vector  sum  of  the  E  and  A  channels  may  have  a  quadrature  component,  as  shown 
in  Figure  4-1. 


2  CHANNEL 
A  CHANNEL 


DIRECT 

MULTIPATH 

TOTAL 

1 

Figure  4-1.  Vector  relationships. 


The  phase  sensitive  detector  used  in  some  radars  to  perform  the 
A/ Z  ratio  function  and  preserve  the  phase  sense  has  the  characteristic: 

|a/z|  -  £  cos 6  4.1 

where  6  m  phase  angle  between  Z  and  A  (ideally  zero). 

For  a  given  target  position,  each  signal,  direct  and  reflected,  is 
weighted  by  the  antenna  gain  corresponding  to  its  angle  of  arrival  for  both 
the  Z  and  A  channels.  Gain  patterns  from  actual  antennas  and  hypothetical 
shaped  beams  were  used.  These  weighted  signals  within  each  channel  are 
vectorially  added  to  yield  the  total  signal  in  that  channel.  Then  |a/e|  is 
computed  using  Equation  4.1  and  the  calibration  look-up  table  is  entered  to 
determine  the  measured  angle.  The  measurement  error  is  the  difference  between 
the  measured  angle  and  the  true  angle.  These  error  calculations  are  performed 
for  a  sequence  of  target  positions  along  the  flight  path  of  the  simulated 
firing  mission. 

The  geometric  definitions  are  shown  in  Figure  4-2.  Note  that  the 
boresight  of  the  antenna  is  maintained  on  the  gunnery  target.  The  angle 
measurement  error  is  calculated  for  the  object  which  is  located  at  target 
altitude  plus  miss  distance  altitude.  Setting  miss  distance  equal  to  zero 
results  in  calculating  the  error  corresponding  to  the  target. 

Three  models  were  used  to  generate  the  input  signals  to  the  angle 
measurement  process: 

1.  A  three-dimensional  facet  model  of  the  terrain  at  the  White  Sands 
Missile  Range.  The  facet  model  is  described  in  Appendix  A. 

2.  A  simplified  two-dimensional  specular  reflection  model  of  the  radar 
two-way  interference  paths. 

3.  A  two-dimensional  specular  reflection  model  of  the  beacon  one-way 


Figure  4-2(b).  Multipath  interference  paths. 


The  signals  for  each  of  Che  radar  (two-way)  paChs  shown  in  Figure  4-2 

are: 


Path 

E  Channel 

A  Channel 

Direct/Direct 

K  (T) 

is 

K  Gj.(T)  .  Ga(T) 

Direct/ Indirect 

K  Gj.(T)  «Gj.(M)p 

K  Gz(T).GA(M)p 

Indirect/Direct 

K  Gj.(M)  •Gj.(T)p 

K  Gz(M).GA(T)p 

Indirect/Indirect 

K  G2(M)2  p2 

K  G_(M)»G  (M)  p2 
L  A 

For  beacon  tracking  (one- 

-way),  the  signals  are: 

Path 

E  Channel 

A  Channel 

Direct 

K  Ge(T) 

K  Ga(T) 

Indirect 

K  G^M)  p 

K  Ga(M)  p 

where 

G  (  )  «  Antenna  gain  as  a 

function  of  angle  off 

boresight 

T  -  Target  angle  off  boresight 

M  *  Multipath  interference  angle  off  boresight 

p  *  Reflection  coefficient 

K  ■  Constant  from  the  radar  range  equation 

E  ■  Subscript  denoting  sum  pattern 

A  *  Subscript  denoting  difference  pattern 


The  magnitude  of  the  voltage  reflection  coefficient  is  an  input  parameter 
to  the  computer  programs  which  calculate  the  angular  error.  Another  input  is 


the  phase  angle  of  the  reflection,  relative  to  the  direct  path,  and  represents 
the  sum  of  the  phase  effects  from  path  length  difference  and  a  complex 
reflection  coefficient. 

The  quadrature  component  is  also  computed  to  determine  its  usefulness  as 
a  measurement  quality  indicator  using: 

QUAD  -  10  l°g10{|f|  sin  s/  4.2 

The  presence  of  a  quadrature  component  indicates  multi; ath  interference; 
however,  relative  phase  angles  of  0  or  180  degrees  between  the  direct  and 
interfering  signal  produce  no  quadrature  component  but  can  result  in 
significant  angle  errors.  The  simple  specular  reflection  model  exhibits  deep 
nulls  in  the  quadrature  component;  therefore,  the  quadrature  component  would 
not  be  a  reliable  track  quality  indicator  if  the  surface  is  a  smooth  plane, 
e.g. ,  a  lake.  The  terrain  of  the  test  range,  however,  generally  appears  rough 
at  microwave  frequencies,  thus  giving  rise  to  diffuse  reflections.  The 
returns  from  the  individual  facets  have  different  amplitudes  and  phases.  The 
vector  sum  of  these  returns  does  not  exhibit  deep  nulls.  Therefore,  in  the 
overland  application  of  MIDI,  the  presence  of  a  quadrature  component  may  be  a 
valuable  track  quality  indicator. 


4. I. 1.2  MIDI  Angle  Measurement 

In  the  current  MIDI  radar,  a  microwave  network  converts  the  basic  E  and 

A  channels  into  (4  +  1  —  ]  and  (4  “  3  Alchannels  (see  Figure  4-3).  The 

2  /2  2  /2 

radar  measures  the  phase  angle,  $  between  the  two  channels.  The  of f-boresight 
calibration  procedure  consists  of  constructing  a  look-up  table  of  $  versus 
angle.  For  a  single  input  signal  as  used  in  calibration,  $  is  computed  from: 


where 

A(e) 
Z(e) 
and  e 


$(Z)  2  arctan  (/2  A(e)/E(e))  , 

voltage  gain  of  elevation  difference  pattern  at  angle  e, 
voltage  gain  of  sum  pattern  at  angle  e  , 
elevation  angle. 


4.3 


The  vector  sum  of  all  received  signals  In  the  £  and  A  channels  Is  used  In 
the  MIDI  angle  measurement.  Two  channels  are  formed: 

A*i+]i  4.4 

1  n 

and  B  *  -J  —  *  4.5 

/2 

Both  £  and  A  are  complex  voltages  and  are  not  generally  In  phase.  The  phase 
angle  $  between  A  and  B  Is  calculated  from: 


4  *  -  arccos 


where  B*  ■  complex  conjugate  of  B. 

The  calibration  look-up  table  Is  used  to  transform  $  Into  the  measured 
angle.  The  measurement  error  Is  the  difference  between  the  measured  angle  and 
the  true  angle. 


4.1.2  COMPUTER  PROGRAM  "COMPARE" 

4. 1.2.1  Program  Summary 

The  program  environment  Is  shown  In  Figure  4-4.  The  Input  file  structure 
is  shown  in  Figure  4-5.  This  structure  which  has  been  used  on  several  other 
MIDI  computer  programs. 

The  program  performs  the  following  operations: 

1.  Reads  in  the  antenna  patterns  and  calculates  the  conventional 
monopulse  and  MIDI  calibration  curves. 

2.  Generates  the  direct  target-to-radar  signal  and  specular  reflection 
signals. 

3.  Weights  each  signal  by  the  antenna  gain  corresponding  to  its  arrival 
angle. 

4.  Computes  the  angle  measured  by  both  conventional  monopulse  and  MIDI 
processing  and  the  measurement  errors. 

5.  Calculates  the  quadrature  component. 
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6.  Repeats  steps  2  through  5  for  each  target  range. 

7.  Creates  an  output  file  containing  plots  of  measurement  error  and 
quadrature  component  versus  range. 

The  program  was  written  in  FORTRAN  and  was  installed  on  the  VAX-11/780 
computer  at  the  Georgia  Tech  Research  Facility,  Cobb  County.  The  output  file 
containing  the  graphs  was  labeled  BITMAP. TBX.  Figures  4-6  and  4-7  are 
representative  plots.  Tabular  data  were  also  output  to  files  FOR007  and 
FOR020. 

4. 1.2. 2  Program  Description 

Program  listings  are  included  in  Appendix  C.  ’’COMPARE"  models  the 
one-way  interference.  "COMPARE  2“  is  identical  except  the  two-way  equations 
are  used.  Each  program  is  separated  by  comments  into  the  following  sections: 

Declarations 
Internal  Data 
Input 

Conversions,  etc. 

Monopulse  Calibration  Curve 
MIDI  Calibration  Curve 
Calculation  Do  Loop 

Generation  of  Direct  and  Reflected  Signals 
Monopulse  Measurement 
MIDI  Measurement 
Plot  Results 

The  listing  is  largely  self-explanatory,  but  some  further  comments  are 
appropriate. 

Internal  Data  -  Contains  four  complex  unit  vectors  which  are  used  to 
convert  real  numbers  into  complex  numbers. 


Input  -  The  first  line  contains  two  integers.  Setting  the  first  equal  to 
"1”  causes  graphs  to  be  generated.  A  zero  supresses  the  plot.  Setting  the 
second  to  ”1"  causes  debug  data  to  be  read  to  an  output  file. 


Conversions 


Within  the  program,  most  angles  are  expressed  in 
milliradians  and  signals  are  expressed  in  volts. 

Monopulse  Calibration  Curve  -  At  large  angles  off  boresight,  the 
calculated  |a/e|  ratio  decreases.  Two  values  of  angle  produce  the  same 

|a/e|  ratio.  The  program  artifically  resolves  this  ambiguity  by  altering  the 
A /E  versus  angle  characteristic  to  be  monotonic.  The  program  outputs  a 

warning  message  and  proceeds.  Obviously  in  running  specific  cases,  the  miss 
distances  and  ranges  must  be  chosen  so  that  off  boresight  scoring  is  not 
attempted  in  this  falsely  calibrated  region. 

Calculation  Do  Loop  -  The  antenna  gain  at  each  angle  is  computed  by  a 
call  to  subroutine  LININT  which  performs  a  linear  interpolation  using  the  gain 
versus  angle  look-up  table. 

Internal  subroutine  ROTATE  performs  a  phase  rotation  of  a  complex  number. 

Plot  Results  -  The  subroutine  call  to  GRAPH1  creates  output  file 
BITMAP. TBX  containing  the  plots,  automatically  scales  the  X  and  Y  axes,  and 
adds  the  X  axis  label,  Y  axis  label,  and  title.  The  call  to  GRAPH2  adds  the 

second  plot  on  the  same  axes.  GRAPH1  and  GRAPH2  call  several  subroutines  from 

the  Graphics  Library  of  the  VAX. 

4.1.3  COMPUTER  PROGRAM  "MULTIDATA" 

4. 1.3.1  Program  Summary 

This  program  is  identical  to  COMPARE  with  one  exception.  Instead  of  a 
simple  internal  specular  model,  it  uses  an  additional  input  file  (FOR030) 
containing  signal  voltages  from  each  of  the  ground  facets  computed  by  the 
three-dimensional  multipath  model.  The  program  environment  is  shown  in  Figure 
4-8,  and  the  structure  of  input  file  FOR030  is  shown  in  Figure  4-9. 

The  input  file  contains  the  "direct/direct"  signal  corresponding  to  the 
path  from  the  radar  to  the  target  and  back  to  the  radar.  The  next  input  data 
point  is  the  signal  from  the  specular  reflection  facet.  The  mean  and  standard 
deviation  of  the  return  from  each  facet  computed  for  a  number  of  pulses, 
typically  200,  are  available  within  the  file.  The  square  of  the  standard 
deviation  represents  the  average  power  received  from  that  ground  facet  for  the 
number  of  pulses. 


MULTI DATA"  environment 
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The  diffuse  returns  from  esch  facet  are  computed  by  the  multipath  model 
for  the  following  signal  paths: 

1.  Direct/indirect  -  The  path  from  the  radar  to  the  target  to  the  ground 
and  back  to  the  radar. 

2.  Indirect/direct  -  The  path  from  the  radar  to  the  ground  to  the  target 
and  back  to  the  radar. 

3.  Indirect/indirect  -  The  path  from  the  radar  to  the  ground  to  the 
target  and  back  to  the  ground  and  finally  to  the  radar. 

The  signals  of  (1)  and  (2)  above  are  equal  and  are  added  together  within 
the  multipath  model  program.  Program  MULTIDATA  separates  the  two  components 
so  that  they  may  weighted  by  the  antenna  gains. 

Program  MULTIDATA  uses  the  direct  signal  and  the  standard  deviations  from 
each  facet  to  calculate  the  total  signal  in  both  the  Z  and  &  channels.  Each 
signal  generated  by  the  three-dimensional  multipath  model  is  weighted  by  the 
antenna  gain  corresponding  to  its  angle  of  arrival.  The  weightings  used  are: 


Path 

E  Channel 

A  Channel 

Direct/Direct 

g£(t) 

Gj.(T)  •  GA(T> 

Direct/ Indirect 

GZ(T)  .  GZ(M) 

GZ(T)  .  Ga(M) 

Indirect/ Indirect 

Gz<M)  *  Gz(T) 

Gz(M)  *  ga(t) 

Indirect/ Indirect 

Gs<“> 

Gj.(M)  •  Ga(M> 
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G( 8)  is  the  antenna  gain  as  a  function  of  angle  off  boresight  and 


Z  *  subscript  denoting  sum  pattern 

A  -  subscript  denoting  difference  pattern 

T  -  target  angle  off  boresight,  and 
M  »  multipath  interference  angle  off  boresight. 

The  total  signals  in  each  channel  are  used  in  exactly  the  same  way  as 
described  for  computer  program  COMPARE. 

4. 1.3.2  Program  Description 

A  program  listing  is  included  in  Appendix  C.  The  program  is  separated  by 
comments  into  the  following  sections; 

Declarations 
Internal  Data 

Input  from  Antenna  Pattern  File 
Conversions,  etc. 

Monopulse  Calibration  Curve 
MIDI  Calibration  Curve 
Interative  Calculation  Loop 

Vector  Summation  of  Returns  from  Facets 
Monopulse  Measurement 
MIDI  Measurement 
Plots  Results 

The  listing  for  MULTIDATA  is  very  similar  to  that  for  COMPARE,  and  the 
same  program  description  information  presented  in  Section  4. 1.2.2  is  also 
applicable  to  MULTIDATA.  The  iterative  calculation  loop  of  MULTIDATA  is 
controlled  by  a  flag  following  the  set  of  data  for  a  given  target  position  as 
shown  in  Figure  4-9.  Therefore,  the  program  can  accept  a  variable  number  of 
target  positions  during  a  scenario  and  non-uniform  steps  in  position  between 
data  points. 


4.1.4  SUMMARY  OF  RESULTS 
X-Band 

Elevation  error  versus  target  range  plots  were  produced  for  each  of  the 
cases  shown  in  Tables  4-1,  4-2,  and  4-3.  The  plots  are  included  as  Appendix 
B.  As  the  target  range  decreases,  the  antenna  elevation  angle  increases,  and 
the  gain  in  the  direction  of  the  multipath  interference  decreases.  The  target 
ranges  at  which  the  angle  measurement  errors  decrease  below  a  threshold  value, 
e.g.,  0.5  milliradians,  were  read  from  the  plots  and  are  shown  in  the  tables. 

The  antenna  pattern  used  in  the  calculation  of  the  data  in  Table  4-1  is 
shown  in  Figures  4-10  and  4-11.  The  envelope  of  the  antenna  sidelobe 
structure  was  used  rather  than  the  detailed  sidelobe  structure.  The  sum 
pattern  sidelobes  are  down  30  dB  at  angles  greater  than  3  times  the  beamwidth 
from  boresight.  This  antenna  pattern  approximates  the  pattern  of  the  current 
MIDI  6  foot  diameter  X-band  antenna. 

The  error  data  from  Table  4-1  are  plotted  in  Figure  4-12.  The  errors  due 
to  multipath  interference  are  less  than  0.5  milliradian  (approximately  0.03 
beamwidths)  at  elevation  angles  greater  than  15  milliradians  (0.7 
beamwidths).  At  an  elevation  angle  of  12  milliradians,  the  errors  are 
approximately  1  milliradian. 

The  results  in  Table  4-2  were  computed  using  the  antenna  pattern  shown  in 
Figure  4-13.  The  difference  pattern  is  shifted  upward  in  elevation  0.33 
degree,  and  the  peak  gain  is  reduced  by  2.5  dB.  The  sum  pattern  remains  the 
same  as  the  one  shown  in  Figure  4-11.  Offset  tracking  is  used.  The  target  is 
tracked  using  a  bias  added  into  the  raonopulse  error  signal  so  that  the 
boresight  (point  of  maximum  sum  pattern  gain)  is  pointed  at  the  target.  The 
errors  of  Table  4-2  are  plotted  with  the  corresponding  non-squinted  cases  in 
Figure  4-14.  The  0.33  degree  (6  milliradian)  upward  squint  permits  the 
antenna  boresight  to  be  lowered  by  5  milliradians. 

The  technique  of  offset  tracking  using  a  normal,  non-squinted  pattern  was 
also  investigated.  The  results  are  summarized  in  Figure  14-15.  The  minimum 
elevation  angle  is  comparable  to  that  achieved  with  the  squinted  pattern. 


TABLE  4-1.  1.3  DEGREE  (22  MILLIRADIAN)  BEAMWIDTH  CASES 


RADAR 

ALTITUDE 

m 


RHO 

MAG/ PHASE  (DEG) 


TARGET 

ALTITUDE 

tn 


MAXIMUM  (m)  RANGE  FOR  ERROR  * 
<1  mrad  <0.5  mrad 


CASE 

RDR  ALT 

m 

RHO 

MAG/ PHASE  (DEG) 

TGT.  ALT 

m 

MISS 

m 

MAX  RANGE 

<0.5  mrad 

(m)  FOR  ERROR 

<0 . 2mrad  <0 . 1  mrad 

i,\»  •!  .C  C 

• 

1 

5 

0.7/90 

15 

0 

6000 

4400 

3900 

2 

5 

0.7/90 

20 

0 

8300 

5600 

5000 

3 

5 

0.7/90 

30 

0 

12500 

8700 

7600 

4 

5 

0.7/90 

50 

0 

20400 

14300 

12600 

5 

5 

0.7/90 

60 

0 

24000 

17200 

15000 

:• 


Figure  4-10.  1.3  Degree  Beamwidth  Antenna  Pattern. 
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Performance  Summary 


Pattern,  Squinted 


TARGET 

ALTITUDE 


Performance  Summary  —  Xband,  Squinted  Difference  Pattern 


ALTITUDE 


Performance  Summary  —  Xband  Offset  Tracking 


K,  Band 
A 


Table  4-3  and  Figure  4-16  summarize  that  performance.  The  0.4  degree 
(6  mllllradian)  antenna  beamwldth  corresponds  to  a  6  foot  diameter  antenna. 
At  elevation  angles  above  4  milllradlans  (0.7  beamwldths),  the  errors  are  less 
than  0.1  mllllradian  (0.02  beamwldths).  When  normalized  to  the  antenna 
beamwldth,  these  results  closely  follow  the  X-band  results. 

Other  Results 

The  reflection  coefficients  predicted  by  the  three  dimensional  facet 
model  using  the  terrain  characteristics  of  an  area  of  the  White  Sands  Missile 
Range  are  shown  in  Table  4-4.  The  conclusions  presented  herein  were  based  on 
a  constant  reflection  coefficient  value  of  0.7  for  all  ranges.  The  reflection 
coefficient  determined  from  the  tests  at  the  North  McGregor  Range  described  in 
Section  4.2  is  also  shown. 

The  relative  phase  angles  between  the  direct  signal  and  the  diffuse 
multipath  Interference  signals  from  the  ground  facets  are  distributed 
throughout  the  region  *ir/2.  Therefore,  the  presence  of  a  quadrature  component 
in  the  elevation  difference  channel  is  useful  in  indicating  the  presence  of 
multipath  interference. 

4.2  MULTIPATH  INTERFERENCE  MEASUREMENTS 

As  was  done  for  the  clutter  analysis,  the  present  MIDI  system  was  used  to 
document  the  magnitude  of  multipath  interference  in  a  desert  environment.  A 
metallic  balloon  was  used  to  measure  the  vertical  lobing  structure 
characteristic  of  multipath  interference,  from  which  estimate  an  of  the 
average  forward  scattering  reflection  coefficient  was  derived.  This 
reflection  coefficient  was  compared  to  the  values  used  in  the  theoretical 
models  to  determine  the  validity  of  the  tracking  error  predictions.  The 
following  paragraphs  describes  the  measurement  procedure  and  the  analysis  of 
the  results. 


TARGET  ALTITUDE 


Performance  Summary  Ka  Band 


TABLE  4-4  REFLECTION  COEFFICIENT  COMPARISON 


TEST 

NORTH 


M  /P8  +  pd 


SIMPLE  FACET 

MODEL  MODEL 


MCGREGOR 


4.2.1  MULTIPATH  INTERFERENCE  MEASUREMENT  PROCEDURE 


The  multipath  measurements  were  preceded  by  a  step  calibration  of  the 
receiver  response  using  10  dB  attenuators.  The  procedure  was  identical  to 
that  used  in  the  clutter  measurements.  Once  calibration  was  completed,  the 
multipath  measurements  proceeded  as  follows.  A  1  m^  RCS  aluminum  balloon  was 
inflated  with  helium  and  attached  to  a  string  marked  off  in  2  foot 
increments.  The  balloon  was  positioned  in  an  area  that  afforded  allow  clutter 
background  while  providing  a  clear  line  of  sight  to  the  desert  surface.  Such 
an  area  was  located  South  south  east  of  the  MIDI  site  at  a  range  of  5800  m. 
Once  the  balloon  was  positioned  at  the  site,  the  antenna  was  pointed  slightly 
above  the  horizon.  Once  measurements  began,  the  antenna  position  remained 
fixed.  Thus,  during  data  collection,  the  only  variable  was  the  height  of  the 
balloon.  This  height  was  varied  from  10  feet  to  100  feet  above  the  local 
terrain  in  5  foot  increments.  At  each  position,  several  seconds  of  range 
gated,  sum  channel,  video  data  were  recorded.  Once  the  vertical  probe  had 
veen  completed,  the  balloon  was  raised  to  200  feet,  the  antenna  was 
boresighted  on  it,  and  several  seconds  of  video  data  were  collected.  This 
calibration  measurement  provided  a  sample  of  balloon  backscatter  power  in  the 
absense  of  multipath  Interference.  Once  the  measurement  was  completed,  the 
balloon  was  moved  slightly  in  range  and  azimuth  and  the  entire  process  was 
repeated. 

4.2.2  MULTIPATH  MEASUREMENT  RESULTS 

The  results  of  the  two  multipath  measurements  are  presented  in  graphical 
form  in  Figures  4-17(a)  and  (b).  The  step  calibration  was  used  to  convert  the 
collected  video  data  into  dB  relative  to  the  maximum  step  calibration  level. 
An  attenuation  level  of  0  dB  corresponded  to  70  dB  above  noise  level.  The 
calibrated  data  at  each  balloon  position  were  averaged,  and  the  averaged 
values  were  plotted  as  a  function  of  balloon  position  above  ground  level.  The 
two  figures  exhibit  a  definite  multipath  lobing  structure.  While  the 
structure  is  not  perfectly  uniform,  it  is  strongly  repeatable,  as  can  be  seen 
by  overlaying  the  two  plots.  At  the  top  of  each  figure  is  the  video  sample  of 
the  balloon  in  free  space.  The  multipath  lobing  structure  clearly  oscillates 
around  this  value. 


The  results  of  the  multipath  Interference  measurement  are  presented  In 
tabular  form  In  Table  4-1.  For  each  lobe,  the  peak-to-null  signal  variation 
and  the  lobe  width  are  recorded.  For  data  set  #1,  the  average  peak-to-null 
variation  was  13.5  dB  and  the  lobe  width  was  14  feet.  For  data  set  #2,  the 
lobe  width  was  the  same  and  the  amplitude  variation  was  12.1  dB.  The  average 
peak-to-null  variation  for  both  sets  was  12.8  dB. 

For  monostatic  (two-way)  multipath  Interference,  the  peak-to-null 
variation  Is  a  function  of  the  magnitude  and  phase  of  the  reflection 
coefficient.  At  the  low  grazing  angles  considered  In  the  present  analysis  the 
phase  of  the  reflection  coefficient  is  very  close  to  180°  for  vertical 
polarization.  The  peak-to-null  (PN)  amplitude  variation  is  given  by 

PN  -  l-  +  2P-t-P  ,  4.7 

1  -  2p  +  pZ 

where  p  -  the  magnitude  of  the  reflection  coefficient. 

This  equation  can  be  solved  for  p  by  converting  it  to  the  quadratic 
equation 

(PN  -  l)p2  -  (PN  +  I)2p  +  PN  -  1  *  0.  4.8 

Equation  4.8  and  an  average  peak-to-null  value  of  12.8  dB  which  corresponds  to 
an  amplitude  ratio  of  4.37  were  used  to  calculate  the  reflection  coefficient 
for  the  present  MIDI  environment;  the  result  was  0.35.  This  value  is  less 
than  that  predicted  by  the  diffuse  multipath  model  and  less  than  the  value 
used  in  the  multipath  Interference  tracking  error  predictions.  At  longer 
ranges  and  lower  elevation  angles,  however,  the  reflectivity  will  be  greater 
than  0.35. 


4.3  CORRELATION  DISTANCE  ANALYSIS 


Multipath  Interference  causes  angle  tracking  errors  when  the  MIDI  is  used 
to  track  a  gunnery  target  during  certain  desired  low  altitude  test 
scenarios.  Similarly,  multipath  interference  causes  unacceptable  errors  in 
the  measurement  of  angle  off  boresight  of  bullets  and  missiles  fired  at  the 


TABLE  4-5.  MIDI  MULTIPATH  MEASUREMENT  RESULTS. 
Data  set  //I 


Peak-to-Null 

Lobe  Width 

Lobe 

(dB) 

(feet) 

1 

13 

— 

2 

19 

15 

18 

3 

14 

15 

10 

4 

9 

10 

11 

5 

10 

15 

(  6) 

Results 

average  *  13.5 

average  *  14 

Data  set  #2 

Peak-to-Null 

Lobe  Width 

Lobe 

(dB) 

(feet) 

1 

8 

15 

15 

2 

22 

15 

12 

3 

10 

15 

11 

4 

9 

10 

10 

Results 

average  ■  12.1 

average  «  14 

average  **14 


-  4.25  m 


Total  Results 


average  ■  12.8 


target.  If  the  errors  in  target  angle  and  projectile  angle  are  highly 
correlated,  they  will  tend  to  cancel  in  the  miss  distance  calculation.  The 
analysis  described  herein  provides  some  insights  into  the  size  of  the  region 
about  the  gunnery  target  within  which  the  projectile  and  target  errors  might 
be  considered  as  being  correlated. 

The  relative  phase  between  the  direct  signal  and  the  reflected  signal  is 
directly  related  to  the  geometric  path  length  difference.  If  these  relative 
phase  angles  for  the  gunnery  target  and  the  projectile  are  the  same,  then  the 
resulting  angular  errors  will  be  correlated  and  will  tend  to  cancel  in  the 
miss  distance  calculation.  The  specific  degree  of  error  cancellation  is 
dependent  upon  the  characteristics  of  the  radar,  principally  transmitter 
frequency  and  antenna  patterns.  This  investigation  consists  of  two  analyses: 

1.  Calculating  the  geometric  path  length  differences.  The  results  are 
expressed  in  meters  and  are  Independent  of  radar  characteristics. 

2.  Determining  the  error  correlation  produced  by  the  path  length 
differences  using  the  characteristics  of  the  current  MIDI. 

4.3.1  GEOMETRIC  PATH  LENGTH  DIFFERENCE  CALCULATION 

The  geometric  definitions  used  throughout  this  analysis  are  shown  in 
Figure  4-18.  The  difference  between  the  reflected  path  length  and  the  direct 


path  length  6j  for  the  gunnery  target  is  given  by: 

61-R1+R2-R  4.9 

where  R1  +  R2  ■  ((hT  +  hR)2  +  G2)*^2  4.10 

and  R  -  ((hT  -  hR)2  +  G2)1/2  .  4.11 

Similarly,  u.e  path  length  difference  for  the  projectile  62  is  given  by: 

6.  -  Rl'  +  R2'  -  R’  .  4.12 


Rl*  +  R2*  -  ((hT  +  hR  +  +  (G  +  Mg)4)1'  4.13 

and  R*  -  ((hT  -  hR  -  M^2  +  (G  +  M^2)172.  4.14 

The  change  A  in  path  length  difference  between  the  projectile  and  the  gunnery 
target  is  then: 


A  -  $2  ”  4.15 

The  quantity  a  is  computed  as  a  function  of  vertical  miss  distance  while 
other  parameters  are  held  constant.  A  listing  of  the  computer  program 
CORRELATION-DISTANCE  which  calculates  and  plots  A  versus  miss  distance  is 
included  in  Appendix  C.  The  program  inputs  are  radar  altitude,  gunnery  target 
altitude  and  range,  and  relative  range  gate  position.  Mg. 


4.3.2  ANGLE  ERROR  CALCULATION 


Angle  errors  are  calculated  for  both  the  target  and  projectile  in  the 
same  way  as  previously  described  in  Section  4.1,  except  as  follows.  Instead 
of  using  a  constant  value,  such  as  90  degrees,  for  the  relative 
phase  ^  between  the  direct  signal  and  the  multipath  interference  signal,  a 
value  is  computed  for  each  gunnery  target  position  and  projectile  position 
using: 

IT6-  4'16 

where 

4  -  phase  angle  of  ground  reflection  coefficient,  taken  as  x  ,  and 

P 

A  *  transmitter  wavelength,  meters. 

Figures  4-19,  4-21,  and  4-23  show  the  relative  error  that  result  from 
the  A  values  of  Figures  4-20,  4-22,  and  4-24,  respectively.  A  reflection 
coefficient  of  0.7  was  used  in  each  case.  Figure  4-25  shows  the  results  for 
the  same  geometry  as  Figure  4-24  but  with  a  reflection  coefficient  of  0.25. 
The  large  error  peaks  occur  when  the  relative  phase  between  the  direct  path 
and  the  reflected  signal  is  a  multiple  of  it.  The  error  budgets  of  0.5 

milliradian  for  ranges  of  4000  meters  or  less  and  0.2  milliradian  for  a  range 
of  10,000  meters  are  exceeded  for  vertical  miss  distances  greater  than  one  to 

two  meters.  The  results  presented  here  tend  to  be  somewhat  pessimistic  in 
that  a  diffuse  reflecting  surface  would  not  exhibit  the  sharp  error  peaks. 

In  conclusion,  the  region  about  the  target  within  which  the  path  length 
differences  were  small  compared  to  a  wavelength  was  not  large  enough  to  be 
useful  for  MIDI  gun  scoring  missions,  which  are  typically  performed  at  ranges 
of  4  km  or  less.  At  ranges  of  10  km  or  greater,  the  technique  may  be 
applicable  for  an  X-band  or  lower  frequency  radar.  The  technique  must  be 
verified  experimentally  before  it  can  be  relied  upon  to  reduce  the  miss 
distance  measure!  nt  errors  in  an  operational  system. 
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Figure  4-19.  Relative  error  at  2,000  m. 
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Figure  4-20.  Path  length  difference  at  2,000  m. 


4.4  CALCULATION  OF  CLUTTER-TO-TARGET  AND  MULTIPATH  I NTERFERENCE-TO-TARGET 


RATIOS 

The  clutter-to-target  ratios  and  multipath  interference  ratios  for 
various  MIDI  II  scenarios  were  calculated.  This  analysis  provided  an  insight, 
expressed  In  the  user  coordinate  system  of  target  range  and  target  altitude, 
into  which  interference,  clutter,  or  multipath,  would  limit  the  radar  tracking 
and  scoring  performance.  Also,  the  effects  of  changing  target  altitude  and 
antenna  beamwidth  could  readily  be  evaluated.  The  model  and  the  computer 
program  that  implements  the  model  are  described  herein. 


4.4.1  INTERFERENCE  RATIO  EQUATIONS 

The  MIDI  II  geometry  is  as  defined  previously  in  Figure  4-2.  The  radar 
target  is,  depending  upon  the  mission,  an  actual  gunnery  target  or  a 
projectile  fired  at  the  gunnery  target.  By  proper  choice  of  radar  cross 
section,  calculations  were  made  for  various  projectiles  and  gunnery  targets. 

The  basic  equations  used  for  the  sum  channel  are: 


PRT  ”  K  GT  aT  * 


PRC  *  K  GC  V 


PRM  "  K  GTGM  p  °T’ 


where 

PRT  *  Power  received  from  target, 
pRC  ■  Power  received  from  clutter, 

PRM  “  Power  received  from  multipath, 
p  «*  Magnitude  of  ground  reflection  coefficient, 

G  -  Antenna  gain,  sum  channel, 

o  ■  Radar  cross  section. 

Subscripts  T,  C,  M  denote  target,  clutter,  and  multipath 
respectively. 
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K  -  Other  parameters  of  the  radar  range  equation  which  are 
constants  for  this  purpose. 

The  clutter  cross  section  is  itself  a  function  of  range  and  was  computed 

from: 

(j  «  R  BM  4R  o  ,  4.20 

o 

where 

R  ■  range  to  clutter  patch, 

BU  -  antenna  azimuth  beamwidth, 

AR  -  range  resolution, 

oq  ■  ground  backscatter  coefficient. 


The  quantities  shown  below  are  derived  from  equations  4.17  through  4.20. 


and 


R  BW  AR  o 

o 
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where  all  quantities  are  as  previously  described. 


4.21 

4.22 


The  basic  equations  for  the  difference  (A)  channel  are  similar  to 
Equations  4.17,  4.18,  and  4.19: 
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where  G.  ■  difference  channel  *>"tenna  gain  and  the  other  quantities  are  as 
A 

previously  defined.  The  expressions  shown  below  are  derived  from  Equations 
4.20,  4.23,  4.24  and  4.25. 


86 


\  --  \  “V  *.  V 


PRCA  _  GC  gac  qc 
PRTA  GT  GAT  °T 


gam  p 


Equations  4.21  and  4.26  are  computed  as  a  function  of  range  while  all 
other  parameters  are  held  constant  and  plotted  on  the  same  axes.  A 
representative  output  is  shown  as  Figure  4-26.  Equations  4.22  and  4.27  are 
similarly  computed  and  plotted;  Figure  4-27  is  a  representative  plot. 

4.4.2  COMPUTER  PROGRAM  "INTERFER_TO_SIG" 

4.4. 2.1  Program  Summary 

The  program  environment  is  shown  in  Figure  4-28.  The  input  file  is  of 
the  structure  previously  shown  in  Figure  4-5,  a  structure  which  has  been  used 
on  several  other  MIDI  computer  programs.  The  other  inputs  requested  by  the 
program  via  the  terminal  are  backscatter  coefficient  in  dB,  target  cross 
section  in  dBsm,  azimuth  beamwidth  in  degrees,  range  resolution  in  meters,  and 
a  title  for  the  plot  output. 

The  program  was  written  in  FORTRAN  and  was  installed  on  the  VAX  11/780 
computer  at  the  Georgia  Tech  Research  Facility,  Cobb  County.  The  program 
creates  an  output  file  labled  BITMAP. TBX,  which  is  subsequently  printed  on  the 
line  printer. 


4. 4. 2. 2  Program  Description 

A  listing  is  included  in  Appendix  C. 
comments  into  the  following  sections: 


The  program  is  separated  by 


Declarations 
Internal  data 
Input  from  file 
Input  from  terminal 


V-  V 


*  .  •  .  v.  •  .  *  .  *  - 


V  V 


AD-A148  467 
UNCLASSIFIED 


METHODOLOGY  IMPROVEMENT  PROGRAM  FOR  THE  SECOND 
GENERATION  MISS  DISTANCE  R.  .  (U)  GEORGIA  INST  OF  TECH 
ATLANTA  ENGINEERING  EXPERIMENT  STATION 
G  V  MORRIS  ET  AL.  30  OCT  84  F/G  17/9 


F/G  17/9 


Figure  4-28  Program  "INTERFER-TO-SIG"  environment 


Conversions 
Calculation  do  loop 
Plot  results 


The  listing  is  largely  self-explanatory,  but  some  further  comments  are 
appropriate. 

The  following  sequence  is  repeated  as  the  target  position  is  stepped  from 
the  minimum  to  maximum  range.  The  program  maintains  the  antenna  boresight 
pointed  at  the  target.  The  angles,  referenced  to  antenna  boresight,  to  the 
clutter  patch  and  to  the  specular  reflection  point  •  are  calculated.  The 
antenna  gain  at  each  angle  is  computed  by  a  call  to  subroutine  DQUAD1NT  which 
performs  a  quadratic  interpolation  using  the  gain  versus  angle  look-up 
table.  The  multipath  interference  ratio,  clutter  to  target  signal  ratio,  and 
target  range  are  stored  in  arrays  MULT  TO  TGT,  CLUT  TO  TGT,  and  TGT  RNG, 
respectively. 

The  subroutine  call  to  GRAPH1  creates  an  output  file  containing  a  plot, 
automatically  scales  the  X  and  Y  axes,  and  adds  the  X  axis  label,  Y  axis 
label,  and  title.  The  subroutine  call  to  GRAPH2  plots  the  second  curve  on  the 
same  set  of  axes.  GRAPH 1  and  GRAPH2  call  several  other  subroutines  from  the 
Graphics  Library  of  the  VAX. 
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SECTION  5 

HARDWARE  FEASIBILITY 


This  section  reviews  those  hardware  subsystems  which  are  key  to  the 
performance  of  the  radar  described  in  this  report.  The  key  subsystems 
reviewed  Include  the  transmitter,  the  signal  processor,  and  the  computer.  The 
objective  of  this  review  is  to  demonstrate  that  a  radar  can  be  assembled  using 
major  subsystems  which  are  known  to  be  off-the-shelf  items. 

For  the  transmitter,  the  goal  is  a  single-pulse  signal-to-noise  ratio  of 
20  dB  at  the  maximum  scoring  ranges  for  each  projectile  with  a  standard 
frequency  and  standard  power  transmitter  tube. 

The  signal  processing  data  rate  was  calculated  in  order  to  estimate  the 
size  and  complexity  of  the  signal  processor. 

For  the  computer  requirements,  the  goal  is  to  Identify  at  least  one 
commercial  product  which  meets  the  requirements  of  MIDI  II  and  which  is  well 
supported  logistically  both  in  hardware  and  software. 

The  data  processing  should  retain  the  essential  features  of  the  current 
MIDI.  The  signal  processing  and  data  processing  discussions  of  this  section 
treat  primarily  (1)  the  new  or  modified  requirements  which  result  from  using  a 
Doppler  radar  with  a  large  number  of  range  gates  and  (2)  upgrading  the 
computing  equipment  to  provide  the  final  vector  miss  distance  calculation 
within  two  minutes. 

5.1  TRANSMITTER  POWER  REQUIREMENTS 

Table  5-1  lists  the  terms  of  the  range  equation  for  three  frequencies 
(10,  17,  and  35  GHz)  for  a  3  foot  diameter  antenna  and  a  transmitter  with  a 
peak  power  of  50  kilowatts.  At  a  range  of  1000  meters,  the  calculated  signal- 
to-noise  ratios  are  36.6,  41.1,  and  48.4  dB,  respectively,  for  these  three 
frequencies.  These  calculated  values  include  13  dB  from  the  pulse  compression 
but  do  not  include  the  15  dB  signal-to-noise  ratio  maximum  improvement  due  to 
the  Doppler  processing. 


TABLE  5-1.  MIDI  II  SIGNAL-TO-NOISE  RATIO  (SNR) 
(o  *  0.001  ■  -30  dBsm) 


TABLE  5-2.  MIDI  II  MATRIX  OF  POSSIBILITIES 


Bullet  scoring  signal  to  noise  ratio. 


35  GHz  &  0.32m 
3  foot  antenna 


Missile  scoring  signal  to  noise  ratio 


0. 001m' 


Figure  5-3.  Signal  to  noise  —  10.5  foot  X  band  antenna. 


Figures  5-1,  5-2,  and  5-3  are  extensions  of  the  data  from  Table  5-1  for 
different  frequencies,  target  sizes,  antenna  diameters,  and  target  ranges. 
The  key  point  is  the  ability  of  a  50  kilowatt  peak  power  transmitter  to 

develop  adequate  (i.e.,  20  dB)  single-pulse,  signal-to-noise  ratios  for 
selected  combinations  of  frequency  and  antenna  diameter. 

Table  5-3  shows  frequency,  coherency,  tube  type,  duty  percent  (%),  price, 
pulsewidth,  peak  power,  and  tube  delivery  time.  There  are  no  significant 
differences  in  the  tubes  required  for  the  different  frequencies.  In  general, 
10  GHz  tubes  are  better  understood  and,  therefore,  easier  to  engineer  into  a 
transmitter  than  either  the  17  GHz  or  the  35  GHz  alternative.  A  50  kW,  10  GHz 
transmitter  with  a  pulse  compression  ratio  of  13  to  1  achieves  the  same  range 
performance  as  a  short  pulse  radar  with  a  transmitter  peak  power  of 
approximately  650  kW.  At  X-band,  1  MU  tubes  exist,  and  a  short  pulse 

transmitter  could  be  used.  At  17  GHz  or  35  GHz,  the  required  peak  power  of  a 
tube  to  implement  a  short  pulse  radar  exceeds  500  kW,  and  tubes  are  not 
readily  available. 

5.2  SIGNAL  PROCESSING  DATA  RATE 

Figure  5-4  is  a  representation  of  the  MIDI  II  sampling  space.  The 
sampling  space  is  composed  of  141  contiguous  range  bins,  and  each  bin's  range 
extent  is  16  feet.  The  total  extent  of  the  range  window  is  2256  feet.  The 

141  samples  are  reduced  by  the  pulse  compression  process  to  141  -  13  *  128 

range  bins.  The  total  energy  is  illustrated  as  contained  in  an  angular  space 
bounded  by  AEL  and  AAZ,  the  3  dB  beamwidths  of  the  antenna. 

The  sampling  space  of  Figure  5-4,  a  projectile  velocity  range  of  500  to 
1000  meters  per  second,  and  a  PRF  rate  of  4000  Hertz  were  used  to  develop  the 
following  reference  values: 

Range  Bin  -  16  feet/4.88  meters 

Projectile  Velocity  -  500  to  1000  m/s 

PRF  -  4000  Hz 

Projectile  Time 

Per  Range  Bin  -  0.005  to  0.01  seconds 


TRANSMITTER  TUBE  REVIEW 


AZIMUTH  -  ELEVATION 


Figure  5-^.  MIDI  II  sampling  space 


Number  of  Samples 

Per  Range  Bln 

Per  Projectile  -  20  to  40  samples 

Figure  5-5  contains  a  high  level  analysis  of  the  signal  functional 
flow.  The  signal  flow  of  this  figure  is  composed  of  these  parts: 

1.  A  transmitter  and  receiver  develop  four  bipolar  voltages  of  the  form 

±1/2  ±  j  A(  )//2  . 

2.  The  four  bipolar  signals  are  in-phase  and  quadrature-phase  sampled 

once  for  each  16  foot  range  bin.  This  requires  a  sample  rate  of 

approximately  32  MHz.  A  nine-bit,  analog-to-digital  (A-to-D)  device 
was  chosen  because  of  their  availability  and  the  attendant  54  dB  of 
dynamic  range. 

3.  The  eight  nine-bit  samples  are  stored  in  a  memory  for  32  PRI's.  The 
structure  of  this  memory  is  shown  in  Figure  5-6. 

4.  The  fast  Fourier  transform  is  applied  to  each  of  the  four  complex 

sampled  signals  for  the  32  PRI's.  The  energy  of  each  signal  is 
resolved  into  Doppler  bins  of  width  125  Hertz.  Because  the  maximum 
Doppler  shift  corresponds  to  60,000  Hertz,  the  PRF  of  4000  Hertz 

undersamples  the  Doppler  by  a  factor  of  15.  This  creates  a  15  to  1 
Doppler  foldover  in  the  FFT  output,  but  the  foldover  is  not  a  problem 
as  the  velocity  of  the  projectile  will  be  accurately  known  before  the 
radar  detects  the  projectile.  This  prior  knowledge  can  be  used  to 
resolve  this  Doppler  foldover  (Doppler  ambiguity). 

The  FFT's  required  per  second  if  the  FFT  is  calculated  for  each 
batch  of  32  PRI's  are: 

FFT's  _  .  RANGE.  ,4  COMPLEX  SIGNALS  .  .4000  PRI. 

SEC  U  1  BINS  ;  K  BIN  ;  ^  32  SEC'  5,1 


FUNCTIONAL 

FLOW 


FUNCTION 

DESCRIPTIONS 


TX  &  RX 


Receiver  output  is  four  signals  of  the 

form;  (1)±j./2+J AE1//2;  (2)±jI/2*AEl/v/2- 
(3)±I/2+;AAZ//2;&  <4)±j  £/2+AAZ/>/2 


Memory 

#1 


RANGE  RESOLUTION  =  l6  FEET  *  4.88  METERS 
SAMPLING  FREQUENCY  =  C/2(AR)  =  9-84  X  10° 

/(2)(16) 

SAMPLING  FREQUENCY  =  32-2  MHz 
EIGHT  A-TO-D* S  REQUIRED  (i&Q  FOR  EACH  OF 
FOUR  SIGNALS) 

BITS  PER  A-TO-D  =  9 

BITS  PER  RANGE  BIN  =  (9)  (8)  =  72 

BITS  PER  PRI  =  (72)(l4l)  =  10,152 

HITS  PER  16  FEES?-  20  @l,00Cta/s  &  40  @50Cto/s 

A-TO-D  DYNAMIC  RANGE  *  54dB 


SIZE  =  (9)(l4l)(8)(32)  =  324,864  BITS 
SAMPLES  PER  FFT  =  32 
OBSERVATION  TIME=  8ms 


SAMPLES  PEIR  FFT  =  32 
DOPPLER  RESOLUTION  =  125  Hz 
RANGE  RATE  RESOLUTION  AT  X-RAND  =  2ms 
DOPPLER  FOLDOVER  <  15  TO  1  (i.e.  60,000/ 
4,000) 


Memory 

n 


SIZE  =  (9) (141) (8) (32)  =  324,864 


Figure  5-5.  Signal  functional  flow. 
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DESCRIPTION 


PULSE 

COMPRESSION 


NON  COMPRESSED  PULSEHXDTH  *»  4l6ms 
COMPRESSED  PULSEWIDIH  «  32ns 
PULSE  COMPRESSION  RATIO  -  13 
PULSE  COMPRESSION  RATION  =  224B 
BIT  RANGE  INCREASE  -  9  +  4  =  13  BITS 


SIZE  -  (13) (128) (8) (32)  =  425,984  bits 


AAZ  -  2  t2n'}(/2  AZ/E) 
AE1  -  2  t2n~  El/E) 


SIZE  *  (13)  (126)  (3)  (32)  »  159,744  BITS 


THRESHOLD  FUNCTION  -  NEAREST  NEIGHBOR 
IF  SRiD  iK{(SR.  j  +  (SR_lfD>  + 

<SR-1,I»I>  +  <SR,D-1>  +  <°>  + 
^R.Dfl5  +  (SR+I,D-I)  +  (SR+I,D) 
+  (SR+l,Dfl)J 


THRESHOLD 


W 


MEMORY  SIZE  =  (9  BITS/RANGE  CELL)  (l4l  RANGE  CELLS/MEAS 
(8  MEASUREWENTS/PRI)  (32  PRl) 

=  324,864  BITS 


Figure  5-6.  Memory  number  1  structure 


TIT  “  70’500  * 

The  Westinghouse,  Inc.,  Very  High  Speed  Integrated  Circuit  (VHSIC) 
contract  calls  for  a  Complex  Arithmetic  Vector  Processor  (CAVP)  which 
computes  a  1024-point  FFT  in  131  microseconds.  A  recent  proposal  by 
another  manufacturer  contained  a  hardware  description  for  a  1024- 
point  FFT  in  50  microseconds. 

The  computational  complexity  controls  the  time  execution  of  an 
FFT.  The  computational  complexity  for  an  N-point  FFT  is  N  log2  N.  A 
1024-Point  FFT  has  a  computational  complexity  of  10,240  while  a  32- 
point  FFT  has  computational  compelxity  of  160.  The  complexity  ratio 
these  two  different  length  FFT's  is  64  to  1.  This  ratio  gives  a  32- 
point  FFT  performance  time  of  approximately  2  microseconds.  The  CAVP 
is  expected  to  produce  500,000  FFT's  per  second.  Also,  this 
computation  rate  is  within  the  capability  of  some  commercial  array 
processors  such  as  the  one  described  later  inthis  section.  The 
components  required  to  meet  the  MIDI  II  FFT  requirements  will  be 
available. 

6.  The  FFT  data  are  stored  in  an  amplitude  -  range  -  Doppler  map  (i.e., 
memory  #2). 

7.  A  pulse  compression  ratio  of  13  to  1  is  applied  to  the  data  of  the 
amplitude  -  range  -  Doppler  map.  A  thirteen-bit  correlation  process 
is  one  method  of  performing  the  pulse  compression.  This  process  has 
been  performed  at  rates  of  100  MHz  for  codes  of  64  bits  by  Georgia 
Tech.  The  number  of  correlations  which  can  be  performed  during  a 
time  period  corresponding  to  32  PRI  is  80,000.  The  required  number 
of  correlations  is  4512  per  32  PRI. 

8.  The  compressed  pulse  data  are  stored  in  an  amplitude  -  range  - 
Doppler  map.  This  map  is  contained  in  memory  #3  whose  structure  is 
illustrated  in  Figure  5-7. 


EDGE  VIEW 


Figure  5-7 •  Memory  number  3  structure . 


9.  The  sum  and  angle  error  terms  are  formed.  The  required  computation 
rate  Is  estimated  as: 

(20  Ops)  (128  Range  Bins)  ,32  Doppler  BinSs 

Req.  Op.  Rate  -  - = - RanS-e  Bln- 

n  r  8  ms 


5.2 


Req.  Op.  Rate  -  10.24  x  106  Ops/sec. 


10.  The  sum  and  angle  error  voltage  are  stored  in  memory  #4. 


11.  A  nearest  neighbor  or  “thumb  tack"  threshold  function  Is  applied  to 
detect  targets.  This  requires  less  than  ten  operations  per  range  - 
Doppler  bin.  There  are  4096  range  -  Doppler  bins  so  the  required 
operations  are  40960  or  less.  The  required  operations  rate  is: 

Req.  Op.  Rate  ■  — — «  ■  512,000  Ops  per  sec.  5.3 

8  x  10 

12.  Not  shown  in  Figure  5-2  are  the  range  trackers.  This  function 
divides  into  the  three  functions  of: 

(a)  Real  time  track  of  the  projectile  target. 

(b)  Near  real  time  sorting  of  the  data  into  subsets  which  are 
associated  with  a  particular  projectile. 

(c)  Near  real  time  fitting  of  a  projectile  trajectory  to  each  subset 
of  projectile  data. 

The  real  time  track  of  the  target  can  be  implemented  in  a  early-late  gate 
tracker  followed  by  double  integration  of  the  error  from  the  early-late  gate 
tracker.  The  tracker  will  use  the  data  from  the  threshold  function  as  input 
to  a  digital  tracker  implemented  in  software.  The  target  tracker  will  remain 
on  the  target  because  of  prior  track  data,  because  of  Doppler  difference 
between  the  projectile  and  the  target,  and  because  of  the  amplitude  difference 
between  the  projectile  and  the  target  returns. 

The  near  real  time  sorting  of  the  projectile  data  will  be  based  primarily 
on  the  projectile  Doppler  shifts.  The  initial  range  separation  between 
projectiles  is  equal  to  or  greater  than  approximately  55  feet  if  the  muzzle 
velocity  is  equal  to  or  less  than  1000  meters  per  second  and  the  firing  rate 


is  equal  to  or  less  than  3600  rounds  per  minute.  For  a  five  second  or  less 
projectile  time-of-f light,  the  differential  velocity  required  for  two 
projectiles  to  pass  each  other  is  11  feet  per  second  or  greater.  The  Doppler 
resolution  of  the  FFT  is  125  Hz,  or  6.8  feet  per  second,  at  a  radar  center 
frequency  of  10  GHz.  Therefore,  the  FFT  output  combined  with  the  firing  time 
of  the  projectiles  will  allow  the  data  set  to  be  subdivided  by  projectile. 
The  near  real  time  trajectory  fit  could  be  executed  as  a  Kalman  Filter.  The 
maximum  data  points  per  projectile  will  be  approximately  150  for  a  2048  foot 
range  data  structure  and  a  projectile  with  a  velocity  of  500  meters  per 
second.  The  trajectory  fit  will  use  150  or  less  data  points  composed 
of  T,  R,  R,  Ag,  and  (Time,  Range,  Range  rate,  Azimuth,  and  Elevation). 

5.3  COMPUTER  SELECTION 

Several  hardware  sets  exist  that  can  perform  the  high  speed  signal 
processing  required  to  reduce  the  front  end  data  to  a  data-stream  which  is 
within  the  performance  envelope  of  a  general  purpose  embedded  computer.  One 
choice  is  the  SPS-1000  array  processor  shown  in  Figure  5-8.  This  machine  is 
available  off-the-shelf  from  Signal  Processing,  Inc.  of  Waltham, 
Massachusetts,  and  it  possesses  the  capabilities  needed  to  do  the  MIDI  II 
front  end  processing.  With  the  addition  of  the  eight,  9-bit,  32  MHz,  analog- 
to-digital  converters,  the  SPS-1000  array  processor  could  perform  the  signal 
processing  functions  at  the  required  data  rates  (i.e.,  32-point  FFT  in  10 
microseconds)  for  the  signal  processing  flow  developed  in  Figure  5-2. 

The  General  Purpose  Processor  attached  to  one  of  the  SPS-1000  ports  could 
be  a  PDP-11  class  machine  with  a  software  development  environment  for  both  the 
PDP-11  and  the  SPS-1000.  This  approach  will  ensure  that  the  computer  software 
can  be  supported  in  the  field  for  both  machines. 

A  General  Purpose  Mini-Computer  of  the  PDP-11  class  is  needed  to  perform 
the  functions  required.  The  three  most  critical  functions  are: 

1.  To  supply  a  well  supported  software  development  environment  for  the 
signal  processing  hardware. 


I  VIDEO 
CHANNEL  #1 


Q  VIDEO 
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Figure  5-8.  SPS-IOOO  processing  architecture 
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2.  To  route  the  reduced  raw  data  from  the  SPS-1000  class  machine  to  a 


disk  or  tape  storage  machine. 

3.  To  perform  the  miss  distance  calculation. 

The  first  function  is  ensured  by  the  wide  range  of  applications  for  which  PDP- 
11  machine  have  been  used. 

Some  characteristics  of  the  PDP-11/44M  are  shown  in  Table  5-4. 


TABLE  5-4.  PDP-11/44M  CHARACTERISTICS 


Instruction 

Move 

Add  1.01  microsecond 

Multiply 

Divide 

Floating  Point  Processor 
Add-32  Bit-Single  Precision 
Multiply-32  Bit-Single  Precision 

Memory 

Cycle  Time-Core 
Read  Access 

Input/Output 
I/O  Rate 


Timing 

1.19  microsecond 

6.38  microsecond 
10.76  microsecond 

8.9  microsecond 
16.2  microsecond 

900  microsecond 
452  microsecond 

<  1.1  Mega-Words  per  second 


The  calculated  I/O  Rate  for  the  MIDI  II  is  15,000  16-bit  words  per 
second.  This  is  well  within  the  1.1  MW/s  of  the  PDP-11/44M. 


The  processing  rates  of  the  PDP-11/44M  are  adequate  to  complete  the  miss 
distance  calculation  for  40  targets  within  a  few  minutes  after  the  test  is 
completed. 

In  Figure  5-9,  a  conceptual  organization  of  the  remote  radar  is  shown. 
This  figure  is  based  on  the  assumption  that  all  the  digitized  data  would  be 
transmitted  via  a  fiber  optics  communications  system  to  the  signal  processor 
van.  Further  analysis  of  the  remote  radar  concept  is  contained  in  Section  6. 

In  Figure  5-10,  a  conceptual  organization  of  the  radar  and  the  signal 
processing  are  shown.  The  signal  processing  concept  selection  was  chosen  to 
meet  these  objectives: 

1.  Transmit  the  full  raw  digitized  data  set  to  the  SPS-100G. 

2.  Process  the  raw  data  set  to  produce  raw  detected-target  data  set. 

3.  Record  the  full  detected-target  raw  data  set  on  the  disk. 

4.  Post  operation,  record  the  detected-target  raw  data  and  the  processed 
data  on  the  tape. 

5.  Select  a  computer  and  software  development  environment  which  allows 
both  development  and  maintenance  of  the  software  for  both  the  PDP-11 
and  the  SPS-1000. 

6.  Make  all  operational  programs  available  in  either  the  memory  of  the 
SPS-1000  or  the  PDP-11. 

7.  Make  quick  load  programs  available  on  the  disk. 


Figure  5-9.  Remote  radar  concept. 


Figure  5-10.  Radar 


SECTION  6 

REMOTE  ANTENNA  CONCEPT 


Without  Doppler  processing,  the  MIDI  II  would  be  limited  to  tracking 
targets  at  3.5  or  more  beamwldths  above  the  ground.  With  Doppler  processing, 
the  MIDI  II  would  be  limited  to  tracking  targets  at  0.8  or  more  beamwldths 
above  the  ground.  Even  with  Doppler  processing,  the  desired  coverage  is  not 
completely  provided  when  the  MIDI  II  is  located  behind  the  gun  or  missile 
launcher. 

A  low  risk  technical  alternative  to  a  multifrequency  radar  for  increasing 
the  coverage  is  the  remote  location  of  the  radar  as  initially  shown  in  Figure 
5-9. 

The  prefered  solution  is  composed  of  these  elements: 

1.  A  coherent  X-band  transmitter. 

2.  A  general  purpose  computer  for  local  and  remote  control  of  the  remote 
radar. 

3.  A  fiber  optics  communication  system  to  route  the  data  between  the 
signal  processing  van  and  the  remote  radar. 


These  elements  are  shown  in  the  concept  drawing  Figure  5-10. 

In  Figures  6-1  and  6-2  a  remote  X-band  system  with  beamwldths  of  2.5  and 
1.25  degrees  is  shown  overlaid  on  the  missile  and  gun  scoring  envelopes  of  the 
MIDI  II  requirements.  The  remote  antenna  is  located  at  11.5  km  for  the 
missile  case  and  2.5  km  for  the  gun  case.  The  coverage  percentage  for  the 
missile  case  is  approximately  50  percent  per  site  location  for  the  large 
antenna  at  one  beamwidth  ground  clearance.  The  coverage  percentage  for  the 
gun  case  is  approximately  90  percent  per  site  location  for  the  large  antenna 
at  one  beamwidth  ground  clearance.  Army  test  personnel  can  easily  construct 
other  examples  for  specific  test  scenarios. 
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SECTION  7 
CONCLUSIONS 


The  major  conclusion  of  this  study  is  that  the  MIDI  I  limitation  when 
operating  against  targets  500  feet  (150  meters)  or  more  above  ground  level  can 
be  improved  down  to  15  meters  for  gun  scoring  and  60  meters  for  long  range 
missile  scoring  by  applying  the  combination  of  radar  techniques  listed  below. 

1.  Pulsed  Doppler  Qtoving  Target  _Indicator)  is  required  to  operate  at 
elevation  angles  less  than  3.5  radar  antenna  beamwidths  above  the 
ground  to  attenuate  the  ground  clutter.  MTI  permits  operation  down 
to  0.8  beamwidths  above  the  ground. 

2.  Errors  due  to  multipath  interference  will  be  less  than  0.03  beamwidth 
when  operating  at  elevation  angles  of  0.7  beamwidth  or  greater  using 
conventional  tracking  and  off-boresight  scoring.  If  an  elevation 
tracking  bias  is  inserted  to  lift  the  antenna  beam  further  off  the 
ground,  operation  down  to  0.5  beamwidths  should  be  possible. 

3.  The  radar  must  have  multiple,  selectable  beamwidths  to  meet  both  the 
low  altitude  operation  and  the  scoring  volume  requirements.  Multiple 
transmitter  frequencies,  multiple  antenna  sizes,  or  a  combination  of 
both  can  be  used. 

4.  A  remote  antenna/receiver/transmitter  should  be  considered  as  a 
method  for  providing  a  favorable  elevation  ge  metry.  A  second 
benefit  is  an  increase  in  signal-to-noise  ratio. 

5.  All  of  the  candidate  systems  can  be  implemented  with  current 
technology. 
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THREE  DIMENSIONAL  MULTIPATH 
INTERFERENCE  MODEL 


THREE  DIMENSIONAL  MULTIPATH  INTERFERENCE  MODEL 


The  Multipath  Interference  Model  la  used  to  calculate  the  signals 
arriving  at  the  radar  via  each  of  the  principal  ground  reflection  paths.  The 
model  software  was  developed  by  John  Peifer  and  Mike  West  in  March  through 
June,  1982,  and  was  revised  by  Mike  West  in  May,  1984  and  Mark  Clinard  in  June 
and  July  of  1984.  The  multipath  returns  parameter  are  calculated  and  output 
for  an  isotropic,  unity  gain  antenna  and  a  target  with  a  unity  cross 
section.  Other  computer  programs  which  use  these  outputs  must  account  for 
radar  system  parameters  such  as  varying  antenna  directional  patterns, 
beamwldths,  etc.  This  appendix  provides  the  following  information: 

1.  A  description  of  the  output  of  the  multipath  model. 

2.  A  description  of  the  information  needed  to  use  the  multipath  model. 

3.  A  synopsis  of  the  mathematical  and  computational  techniques  used  to 
produce  the  output  parameters. 

The  multipath  model  considers  a  terrain  composed  of  24  x  175  square 
facets  of  equal  area  with  sides  of  88  meters.  Indexing  the  facets  with 
integers  specifying  the  latitudinal  and  longitudinal  coordinates,  the  model 
references  an  input  file  to  provide  parameters  characteristic  of  each 
individual  facet  (terrain  height,  terrain  type,  etc.).  Figure  A-l  shows  the 
faceting  technique  for  a  similar  40  x  40  facet  terrain.  Figure  A-lc  shows  how 
the  radar  site,  a  target,  and  the  surface  in  between  the  underlying  facets 
would  typically  be  located  with  respect  to  the  faceted  terrain. 

The  most  fundamental  physical  parameter  describing  the  reflection 
properties  for  the  terrain  facet  is  its  dielectric  constant.  Given  the  user's 
information  specifying  a  wet  or  dry  scenario,  the  model  considers  six 
different  types  of  terrain  that  may  characterize  each  facet: 


Al.  M.  S.  West,  "A  Multipath-Clutter  Model  for  SAMS,”  Final  Report  on  Contract 
F33615-81-C-1487,  Georgia  Institute  of  Technology,  Engineering  Experiment 
Station,  December  1983. 


1.  Sea 

2.  Soil 

3.  Grass 


4.  Crops 

5.  Trees 

6 .  Sand 


The  dielectric  constants  were  experimentally  determined  for  the  X-band 
frequency  range.  They  are  shown  in  Table  A-l.  The  Fresnel  reflection 
coefficient  pQ  is  evaluated  for  either  horizontal  or  vertical  polarizations: 
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simp  -  (e^  -  cos  ip) 
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simp  +  (er  -  cos  ip) 


for  vertical  polarization 


for  horizontal  polarization 


where: 


the  incident  angle  on  the  facet 
the  dielectric  constant 


With  the  scenario  geometry,  the  terrain  characteristics,  the  operating 
frequency,  the  transmitted  power,  the  gain  of  the  transmitter,  and  target 
established,  the  multipath  program  computes  the  following  output: 


The  direct  return  from  the  target:  (real  number). 

The  range  and  elevation  angle  from  the  receiver  to  the  specular  facet 
at  which  the  reflected  RF  energy  from  the  target  to  the  surface  to 
the  receiver  obeys  Snell's  Law. 

The  returns  from  the  target  to  the  receiver  via  two  different  paths, 
both  of  which  involve  the  specular  facet:  (complex  number). 


Note  that  all  returns  are  in  terms  of  voltages;  therefore,  the  returns  are 


Note  tnat  an  returns  are  m  terms  ot 
proportional  to  1/R^  instead  of  1/R  . 


TABLE  A-l  DIELECTRIC  CONSTANTS 


a.  The  direct-indirect  path  composed  of  the  ray  from  the  radar  to 
the  target,  the  ray  from  the  target  to  the  specular  facet,  and 
the  returning  ray  to  the  receiver.^ 

b.  The  indirect-indirect  path  composed  of  the  transmitted  ray  from 
the  radar  to  the  specular  facet,  the  reflected  ray  to  the  target, 
and  the  identical  return  path  taken  by  the  direct-indirect  path 
back  through  the  specular  facet. 

4.  The  diffuse  multipath  return  coefficients  from  the  individual  facets 
for  the  indirect-indirect  and  direct-indirect  paths  provided  that  the 
sub-facets  satisfy  the  following  criteria  enumerated  by  West: 

a.  The  sub-facet  is  sufficiently  illuminated. 

b.  The  radar  echo  of  the  sub-facet  is  not  clearly  separated  from  the 
radar  echo  of  the  target. 

c.  The  facet  is  not  shadowed  by  surrounding  facets  as  illustrated  in 
FigureA-^. 

5.  The  range,  elevation  angle,  and  azimuth  angle  to  each  sub-facet  from 
the  radar  site  from  which  the  diffuse  multipath  is  calculated. 

Item  #4,  the  diffuse  return,  is  given  in  the  form  of  the  mean  of  the  Gaussian 
phase  distribution  of  the  return  voltages  and  the  standard  deviation  of  the 
distribution;  thus,  the  standard  deviation  for  the  two  paths  is  a  measure  of 
the  returned  voltage.  (Both  the  mean  and  standard  deviation  are  complex 
quantities. ) 

Figure  A-3  illustrates  the  three  possible  return  paths. 

Several  inputs  which  are  qualitatively  described  here  are  needed  to  use 
the  multipath  model.  First,  a  file  describing  the  terrain  and  each  of  its 
facets  must  be  accessed  to  provide: 

1.  The  longitudinal  and  latitudinal  dimensions  of  the  terrain  in  units 
of  facets  along  with  scaling  factors  and  offsets. 


t* 


2.  The  three  dimensional  normal  vector  the  surface  of  each  facet. 

3.  The  terrain  height  of  each  facet. 

4.  The  terrain  type  of  each  facet. 


A  user  input  file  provides  three  vectors:  the  position  of  the  target,  the 
velocity  of  the  target,  and  the  velocity  of  the  receiver  (normally  set  equal 
to  [0,0,0]). 

Each  time  the  multipath  program  is  run  for  a  new  scenario,  the  user  is 
prompted  for  these  mutable  conditions: 

1.  Horizontal  or  vertical  polarization? 

2.  Wet  or  dry  terrain? 

3.  The  receiver  half-bandwidth? 

4.  Site  height  is  meters  above  the  terrain  surface? 

5.  RMS  surface  deviation  for  the  entire  terrain? 

The  location  of  the  radar  site  is  currently  set  in  the  middle  of  the  24-facet- 
wide  field  at  1000  meters  latitude  and  276  meters  down-range.  The  radar  site 
height  is  set  by  the  coinciding  facet  height  plus  the  input  site  height,  as 
described  in  item  #4  above. 

The  qualitative  outputs  of  the  multipath  model  have  been  described  along 
with  the  inputs  needed  to  use  the  program.  The  exact  equations  used,  the 
exact  criteria  for  evaluating  contributions  to  the  multipath  return,  and  the 
particular  parameters  and  results  of  the  MIDI  analysis  are  then  ready  to  be 
specified. 

The  direct  return  coefficient  is  calculated  from  the  radar  equation: 
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where : 


kj  *  the  Doppler  shifted  wave  number  at  the  receiver, 

R  -  the  line-of-site  range  from  the  receiver  to  the  target, 
Pt  ■  the  transmitted  power, 

Gt  -  the  gain  of  the  transmitter  antenna. 
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In  the  runs  for  the  MIDI  analysis,  the  operating  frequency  was  10.125  GHz. 
Typically,  the  target  was  approaching  the  receiver  at  a  velocity  of  100  m/s. 
The  transmitter  antenna  gain,  an  internal  variable  of  the  program,  was  set 
equal  to  unity. 

For  the  return  from  the  specular  facet,  the  model  first  locates  the 
specular  facet  at  which  the  RF  energy  is  reflected  according  to  Snell's  Law  by 
considering  the  coordinates  of  the  receiver  and  the  target.  The  line-of-sight 
range  from  the  receiver  to  the  specular  facet  and  the  elevation  angle  at  the 
receiver  down  to  the  specular  facet  is  output.  (The  elevation  angle  is  equal 
to  the  complement  of  the  off-normal  incident  angle  and  the  reflected  angle.) 
The  direct-indirect  and  the  indirect-indirect  return  coefficients  are  then 
computed. 

The  direct-indirect  and  indirect-direct  returns  are  equivalent  for  the 
case  of  an  isotropic  antenna  as  assumed  in  the  MIDI  analysis  and  are  summed 
together  for  the  output.  The  summed  direct-indirect  (DI)  coefficient  is 
calculated  from  the  radar  equation  and  the  reflection  coefficients  described 
below: 

DI  -  (2pv  pQ  prms  KGt/R2)exp(jkAR)  (A-4) 

where: 

<  *  X/Pt/(4ir)3 

k  ■  the  transmitted  wave  number  (w/c) 

AR  ■  (receiver  to  the  specular  facet  range)  +  (target  to  the  specular 
facet  range)  -  R 

Pv  -  the  vegetation  reflection  coefficient  (0  <  pv  <  1  and  »  1.0 
for  sea,  soil,  and  sand.) 

pQ  *  the  Fresnel  reflection  coefficient 

prms  *  exp  (“2(koh  8ln 

a.  *  the  user-input  rms  surface  deviatior.. 


All  of  the  terrain  used  for  the  MIDI  analysis  was  of  the  type  "sand." 
The  input  rms  surface  height  deviation  was  0.2  meters.  The  terrain  was  dry 
and  the  radar  transmitted  vertically  polarized  energy.  The  indirect-indirect 
(II)  specular  return  coefficient  is  calculated  from: 


II  -  <PvP0Prm8)2  <(Gt/R2)exp(jkAR)  . 


(A-5) 


The  consideration  of  the  diffuse  multipath  return  involves  a  much  more 
complex  analysis  of  the  geometry  of  the  model.  All  of  the  facets  inside  a 
rectangle  whose  sides  are  defined  by  the  edges  of  the  terrain  and  the 
downrange  coordinates  of  the  radar  site  and  the  target  are  initially 
considered  for  possible  contributions  to  the  multipath  return  at  the 
receiver.  First,  criteria  are  applied  to  decide  which  of  the  88  x  88  square 
meter  facets  will  be  subdivided  into  smaller  subfacets  so  as  to  obtain  better 
approximations  for  Doppler  shifted  returns.  Secondly,  criteria  are  applied  to 
each  facet  or  subfacet  to  determine  if  its  contribution  to  the  diffuse  return 
is  negligible  or  is  to  be  evaluated.  Finally,  the  location  of  each  evaluated 
subfacet  or  facet,  its  DI  return,  and  its  II  return  are  output  in  the  form  of 
range,  elevation  angle,  azimuth  angle,  and  the  standard  deviation  and  mean  of 
the  Gaussian  phase  distribution. 

Ideally,  Figure  A-4  shows  the  facets  which  will  be  considered  for  diffuse 
multipath  evaluation,  but  the  program  model  used  here  considers  (and  perhaps 
tosses  out)  a  much  wider  area  for  the  prospective  facets  as  described 
above. Still,  the  facets  darkened  in  Figure  A-4  for  the  general  case  of 
a  radar  site  and  target  placed  in  a  faceted  terrain  will  be  the  primary 
contributors  of  diffuse  multipath  return. 

The  model  proceeds  to  loop  over  all  the  facets  in  the  rectangular  area 
between  the  site  and  the  target,  evaluating  them  one  at  a  time.  Each  facet  is 
divided  into  25  subfacets  of  equal  area  to  be  evaluated  in  series  along  with 


A2.  S.  P.  Stuk,  "TAC-ZINGER  Clutter  and  Multipath  Models  for  Track-While-Scan 
Systems,"  Software  Documentation  on  Contract  F49620-79-C-0222,  Georgia 
Institute  of  Technology,  Engineering  Experiment  Station,  October  1980. 
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the  remaining  facets  if  its  center  is  located  at  a  line-of-site  range  greater 
than  3000  meters  from  either  the  target  or  the  radar  site.  Thus,  if  large 
Doppler  frequency  variations  occur  over  the  area  of  the  facet,  it  is  to  be 
subdivided. 

As  the  facets  and  subfacets  are  evaluated  serially,  their  contributions 
are  neglected  if: 

1.  The  facet  or  subfacet  lies  outside  of  the  cone  defined  by  the 

receiver's  half-beamwidth. 

2.  The  line-of-sight  (LOS)  receiver-to-target  range  divided  by  the 

product  of  the  LOS  receiver-to-facet  range  and  the  LOS  target-to-facet 

— Q 

range  is  less  than  1  X  10  . 

2  2 

3.  tan  g/tan  3q  <  50 
where : 

3  -  the  angle  between  the  vector  normal  to  the  facet  surface  and  the 

bisector  of  the  angle  between  the  LOS  vectors  from  the  facet  to  the 
receiver  and  from  the  facet  to  the  target. 

3  ■  0.2  (for  our  model) 

o 

If  a  subfacet  or  facet  is  not  eliminated  by  any  of  these  criteria,  the  model 
proceeds  to  calculate  its  contribution  to  the  diffuse  multipath  return. 

First,  the  diffuse  reflection  coefficients  must  be  calculated  for  each  of 
the  two  paths  DI  and  II.  A  geometric  mean  factor  involving  the  incident 
angles  shown  in  Figure  A-5  as  and  ^  an<*  the  rMS  surface  height  deviations 
is  first  calculated: 

PD  •  (1  -  -  »,22)  • 

0 

Psl  »  exp(-2(kohsini(»1)‘  ) 

Ps2  -  exp(-2(kohsln<|>2,2)  . 
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where: 


(A-6) 


The  diffuse  part  of  the  model  assumes  a  vertical  polarization, 
reflection  coefficient  for  the  D1  path  then  follows: 


X  0oFD2|(‘>o‘'v,2|Ft‘'t  *  SIDE 


where: 


LOS  range  from  the  facet  to  the  receiver 


R2  ■  LOS  range  from  the  facet  to  the  target 

SIDE  *  the  length  of  the  side  of  the  subfacet  or  facet 
2  2  2 

aQ  m  exp  (-tan  g/tan  gQ)/tan  gQ 
Ft  ■  the  antenna  directional  gain  pattern  (set  *  1) 


The  reflection  coefficient  for  the  II  path  is: 

»oFD2|»o,»vl(SIDB)2  Pt  <Ft>2!o„_»wl 

11 - - - 2 - 572 - 1 - 

(RiR2)Z(4ir):)/Z 


The  diffuse 


(A-7) 


(A-8) 


where: 


pQ  *  the  Fresnel  reflection  coefficient  for 

p  *  ■  the  Fresnel  reflection  coefficient  for  . 

o_  y2 


Given  the  voltage  reflection  coefficients  for  the  two  paths  weighted  by 
the  transmitter  pattern  and  the  area  of  the  glistening  subfacet,  the  model  can 
calculate  the  reflected  power  from  the  facet  (or  reflected  voltage). 

The  multipath  model  assumes  that  the  reflected  voltages  are  distributed 
in  phase  in  the  shape  of  a  Gaussian  curve.  The  model  uses  a  Gaussian  weighted 
random  number  generator  to  distribute  the  phases  of  200  samples  of  the  return 
from  the  subfacet  being  evaluated.  Statistics  a’-*?  then  extracted  from  the  set 
of  samples  to  provide  the  mean  of  the  Gaussian  phase  distribution  and  the 
standard  deviation.  Thus,  the  standard  deviation  of  the  Gaussian  is  a 
relative  measure  of  the  energy  inside  the  envelope  and  reflected  from  the 
subfacet.  The  LOS  range  to  the  center  of  the  subfacet,  the  elevation  angle, 
and  the  azimuthal  angle  with  respect  to  the  receiver  are  then  output  along 


with  the  statistics  for  the  DI  and  II  paths.  The  total  number  of  subfacets 
and  facets  evaluated  for  the  diffuse  multipath  calculations  are  also  output 
with  the  return  data. 

Options  are  available  and  were  developed  by  Mark  Clinard  to  map  the 
diffuse  return  contributions  in  decibels  on  a  printout  of  the  24  x  175  array 
representing  the  faceted  terrain.  Along  with  the  total  contribution  of  each 
facet  being  mapped,  a  map  is  provided  detailing  the  number  of  subfacets  that 
have  evaluated  inside  each  facet. 
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This  appendix  contains  the  angular  error  versus  range 
curves  for  the  cases  shown  in  Tables  B-l,  B-2,  and  B-3. 

Most  cases  contain  curves  for  both  conventional  monopulse 
processing  and  MIDI  processing.  Each  case  is  individually 
labeled . 
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f:i 
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L 
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i-:- 
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TABLE  B-2 

1.3  DEGREE  (22  MILLIRADIAN)  BEAMWIDTH  CASES,  0.33  DEGREE  SQUINT 


CASE 

RDR  ALT 

M 

RHO 

MAG/ PHASE  DEG 

TGT.  ALT 

M 

MISS 

M 

MAX  RANGE 

<1  MR 

(M)  FOR  ERROR 

<0.5  MR 

50 

5 

0.7/90 

15 

0 

2000 

<1000 

51 

5 

0.7/90 

15 

-5 

2000 

1600 

52 

5 

0.7/90 

20 

-5 

3000 

1900 

53 

5 

0.7/90 

30 

0 

4000 

3000 

54 

5 

0.7/90 

30 

-5 

4400 

2900 

55 

5 

0.7/90 

35 

-5 

5100 

3300 

56 

5 

0.7/90 

60 

-10 

8900 

5700 

57 

5 

0.7/90 

50 

0 

6900 

4400 

.  V  V  V 

- 

TABLE  B-3 

' 

0.4  DEGREE  (6 

MILLIRADIAN) 

BEAMWIDTH  CASES 

*  *  - -  j  -  •  * *< 

CASE 

RDR  ALT 

RHO 

TGT.  ALT 

MISS  MAX  RANGE 

(M)  FOR  ERROR 

.  •  ..  7 

M 

MAG/ PHASE  DEG 

M 

M  <0.5  MR 

<0.2  <0.1 

ELEVATION  ERROR 


CLClMTtON  cccoe 


QOO  rttT.3,  TOT  ALT-SO.OMO-.^SG  OCG,  l.JOCO  Bu.nlSS-IS 


000.rH.T- 3,  TGT_/M_r-30,  OMO- .7,-3«  OCO,  I  3  0C0  on,  mss-- IS 


MONO 


CASE  9 


ooo.otT.o#TGr.«i.T.se,PHO...»r30  ocs.i.iocb  eo,mss-0  MONO _ ,  CASE  10 
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MONO 


CLOMTIQN  EBOOfl 


CUVMItON  ERROR 


OOO.fM.f.3,  TGT_ALt>t9,QH0>.7sM  OCG,  130C0  BU,nlSS-B 


WJffi 


QOO.ALT.3, 


>kV« 


count  ion  cpbop  -  m  euumoN 


APPENDIX  C 

COMPUTER  PROGRAM  LISTINGS 

1.  "COMPARE" 

2.  "MULTIDATA" 

3.  "CORRELATION-DISTANCE" 
k.  "MCORELATE" 


5.  "INTERFER  TO  SIG" 


non  n  n  o  r>  n  n 


PROGRAM  COMPARE 


COMPARE  C0NVEN1 IQNAL  MONOPULSE  AND  MIDI 

DIMENSION  SUM  DB(40), DEL JOB < 40 ) * ELANQ_DEQ ( 40 ) . 

CAL_MONO<40>. 

ELANO_MR ( 40  > ,  CAL J1IDI ( 40 ) , 

SUM_VOLT (40)# DEL_VOLT  <  40 ) 

DIMENSION  QUAD ( 200 > 

DIMENSION  TOT JW40< 200 >»  ERROR  J10N0< 200).  ERROR  J1IDI( 200) 

CHARACTER*80  TITLE 
CHARACTER -*40  LABX 
CHARACTER -*40  LABY 
CHARACTER«40  LARY1 
CHARACTERS  TYPEJ.INE 

COMPLEX  SUM  SIOJ)IR.  SUM_SIO  KUlT, DEL  SIO  DIR. DEL  SIO  MULT 
COMPLEX  RHOjCPLX,  A.  B.  SUM_TOTAL,  DEL_TOTAL 
COMPLEX  DEL _ROTAT£U,  UNIT_IN,  UNITjQUAD 

REA!  IN_PHASE#  MISSJDIST 

U4IERNAL  DATA 

DATA  SUMjSIOJDIR/<  1.  .0.  )/ 

DATA  DELJBI0J)IR/<1.  .  0.  )/ 

DATA  UNXT_XN/(1.  .0.  )/ 

DATA  UNIT  QUAD/ ( 0.  .  1.  )/ 

LARX-' TARGET  RANOE  —  M' 

LARY- 'ELEVATION  ERROR  —  HR' 

LARY1- 'QUADRATURE  COMPONENT  --  DB ' 

TYPEJLINB-'DASH' 

INPUT 


READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 


1009)  IPLOT.  IBUO 
1006)  (SUMJDB<  I).  1-1,40) 
1006)  <DEL_PB<  I).  1-1,40) 

1006)  <ELANO_DEO< I ),  I-l»  40) 

1007)  PHASE  JWOJDEO 
1007)  RHO 

1007)  TOT_ALT 
1007)  TOT_MAX  RNO 
1007)  TOT_MIN_RNO 
1007)  RNOJBTEP 
1007)  RDR_ALT 
1007)  MISSJJIST 
1100)  TITLE 


CONVERSIONS.  ETC 


8QR  T2*SQRT ( 2.  ) 

PI-4*ATAN(1.  ) 

P I  jDVfcRJB-P  1/2. 

DEC_T0J1R-PI/.  ISO 
SXCN— 1. 

DO  9  1-1.40 

IF(DELJ)BCI).  LT.  0.  )  SION-1. 

SUM_VOLT ( I )  - 1 0.  **  (  SUM _DB  ( I )  /20.  ) 
DEL_VOLT( I )-3I0N*10.  ••(DEL  DB(I)/20.  ) 
ELANC_MR  ( I )  -ELAM©  J)EQ  ( I )  *DEG_TO_MR 
CONlINUE  lfi? 


nnn  nnn  non  no  non  non  non 


MONOPULSE  CALIBRATION  CURVE 
DO  10  1*1*40 

CAL_M0N0 ( I ) -DEL_V0LT < I ) /SUM  VOLT (I ) 

10  CONTINUE 

PRfcVfcNT  CALIBRATION  CURVE  FROM  BEINO  DOUBLE  VALUED 
DO  15  1*1.39 

IF(CAL_J!0N0(  I).  LT.  0.  )  00  TO  15 
IF (CALJ1QMQ(  I+l ) .  OT.  CAL  MQN0(  I  >  >  00  TO  15 
CAL  J10N0( I+l )* CAL  MONO ( I >+.  0001 
URiTE(6> 1001 )  I 

15  CONTINUE 

DO  16  1*40.  2.-1 

IF(CAL_J10N0(I).  OT.  0.  )  00  TO  16 
IF(CAL_KQN0(I-1).  LT.  CAL_hONO(I))  00  TO  16 
CAL_MONO<  1-1 )  -  CAL_f10N0  <  I )  -.  0001 
URiTE(6.  1008)  I 

16  CONTINUE 

IF ( IBUO.  LT.  1)  GO  TO  25 
WR1TEC7.  1012) 

DO  22  1*1.40 

WR1TF<7. 101 1 )  CAL  MONO ( I ) , ELANOJiR ( I ) 

22  CONTINUE 

25  CONTINUE 

MIDI  CALIBRATION  CURVE 

DO  20  1*1.40 

CAL_MJDI  <  I  )“2.  ♦■AT AN <  C AL_MONO <  I  )*SQRT2> 

20  CONTINUE 

IF <  IBUO.  LT.  1)  CO  TO  24 
WRITE(7.  1013) 

DO  23  1*1.40 

WR|TE(7.  1011)  CAL_MIDI ( I ) »  ELANO  MR(I> 

23  CONTINUE 

24  CONTINUE 

PRELIMINARIES 

PH/.SE_ANG_RAO-PHASE_ANQ_PEG*DEO_TO  MR  / 1 000. 

NTOINIS-O 

TCT JRNC ( 1 )— TOTJMIN  RNG-RNC_STLP 

NSTfcPS-INI  (  (TCTJ1AX_RNG-TCTJ1IN_RNQ>/RN0__STEP>  +  1 

TABULAR  DATA  HEADERS 

WRITE  (20.  1002)  TITLE 
URITE  (20.  1003) 

BEOIN  ENORMOUS  DO  LOOP 

DO  30  1-l.NSTEPS 

CENJ-RATION  OF  DIRECT  AND  REFLECTEO  SIGNAL 

TCJ_RNG(I+1)-T0T  RN0(I)+RN0  STtP 
ANC_0FFJB0RESICHT-1000.  *MISS  DIST/TOT  RN0(I+1> 

AN1 _ANC_DIRECT-ANQ_OFF  BORESfCHT 

CALL  LININT (ELANO  MR,  SUM_VOLT,  40. ANT  ANG_DIRECT,  GAIN  SUM, 
1  IVIEWO) 


uuu  uuu  ooo 


CALL  LZNINT(EI  ANO_MR,  DEL  VOLT#  40#  ANT_ANOJDIRECT#  0AINJDEL# 
IVIEWO) 

IF (  IVIEUO.  NE.  0.  )  WRITE <6,  1009)  IVIEWO 

T«T_AN0JDIR-1000.  *<T0T  ALT-RDR_ALT+MISSJDIST ) /T0T_RN©< 1*1 ) 

TOT _AN©_MULT--1000.  *(T0T  ALT4rtDR_ALT+MISSJ>IBT)/T0T_RN0(  1+1 ) 
ANT  _JANQ_MULT-TQT_ANQ_MULT+ANO_OFF_BOREBIQHT-TQT_AN©_DIR 

CAI  L  LININT (ELANO  MR# SUM  VOLT, 40#  ANT_ANO_MULT,  VOLT_MULT, 
IVIEW1 ) 

SUMJBIC J1ULT-V0LT  J1ULT*RH0*UNI T_IN 
CAI  L  ROTATE (SUH  SI0_MULT,  PHASE_ANO JRAD ) 

CAI  L  LININT (EL ANOJMR.  DELJVOLT,  40.  ANT_ANO_MULT,  VOLT_MULT, 
IVIEW2) 

DEL  JBIC  J1ULT-V0LT  J1ULT*RH0#UNIT_IN 
CAI  L  ROT  ATE  (  DEL  _JB  1 0  J1ULT .  PHASE  _ANO_RAD) 

TOTAL  SUM  ANO  DIFFERENCE  SI0N4>  S 

IF(  (  1VIEW1+IVIEW2) .  OT.  1)  WRITEI6,  1080)  IVIEW1.  IVIEW2,  I 
SU  BTOTAL -SUM  8 I C_D I R*OA I N_SUM+SUM_S I Q_MULT 
DEL  _TOT AL“DEL_S  I C_D I R*0A I N_DEL+DEL JB 1 0_MULT 

tlOhlOPULSE  MEASUKEHENT  ANO  CALCULATION  OF  QUADRATURE  COMPONENT 


81»:ABS(SUM_T0TAL  ) 

D l -CARS ( DEL_TOTAL > 

Rl-REAI  ( SUH_TOTAL* CON JO < DEL_TOTAL ) ) 

I M  J»HA8E-R 1/(81 *+2 ) 

IF(  INJ»HASE.  OE.  1.  2)  IN_PHASE«1.  2 
IF(IN_fHASE.  LE. -1.  2)  IN_PHASE— 1.  2 
Ai<!C— ACQS(R1/(D1-*S1 )  ) 

QUAO(I+l)-  D1*SIN(AN0)/S1 
QUAD( 1+1 )-QUAD( I  +  l >**2 
IF(QUAn<I+l).  LT.  .  0001)  QUAD( 1+1 )— .  0001 
QUAD( I+l )-10.  * ALOQ 1 0 ( QUAD ( 1+1 ) ) 

CALL  LININT (CAL_MONO.  ELAN©_MR,  40.  IN_PHASE# 

ANT_T0T_AN0  MES,  IVIEW3) 

TOT  AN0JC8-ANT_T0T  AN0J1ES-AN0  OFF JBORES I OHT 
IFIIBUO.  EQ.  1)  WRITE (7.  1020)  IN  PHASE#  TGT_ANQ__MES 
ERROR _ftONO(  1+1 )  «TQT_ANQ_MES 

NIDI  MEASUREMENT 

D*-L_P0TATED-DEL_T0TAL/3GRT2 
CALL  ROTATE  (DEL _ROTATED,  PI_OVRR  2) 

A— BUH_TOT AL/2.  +DEL  ROTATED 
B-8UM_T0TAL/2.  -DEL  ROTATED 
DUM1*REAL ( A*CONJ© (B )  ) 

DUM2-C ADS ( A ) >CABS ( B ) 

I F  <  DOM 1 .  GE.  DUM2)  00  TO  27 
C— AC03  (  DUN1/DUH2 ) 

IF(C.  CE.  2.  )  C-2. 

IF(C.  LE.  -2.  )  C— 2. 

00  TO  26 

USE  LAST  VALUE  OF  C 

CONTINUE 

UR1TE(6.  1016)  I 

IF ( IOUO.  EG.  1 )  WRITEI7,  1016)  I 

CALL  LININT(CAL_M7DI#  ELANO _MR,  40,  C,  ANT__T0T_AN0_ME8,  IVIEW4) 
TOT  JNNO  J1ES-ANT  TOT  _AN0_J1ES- AN©  JDFFJB  OR  E8 1  OHT 
IF(dVIfcW3+IVIEW4>.  OT.  1)  WRI1E(6, 1081)  IVIEW3#  IVIEW4#  I 
ERROR _ftIDI  (  1  +  1 ) -TCT_ANO_MES 
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.v\;  .*• 


•*  ,  •*.  **.  ^ .*•*-•*.  *\-\ 


IF ( IBUO.  EQ.  1 )  WKlTE<7, 1021 )  C.  TOT _ANQ_MES 

npoinis-nhoints+i 

UR 1 1 E ( 20, 1004)  T0T_RN0(I+1 ),  TOT _ANQ_DIR, 

3  ERR0RJ10N0U+1), 

1  ERROR  MIDI(I  +  1>,QUA0U  +  1),  IVIEU1, 

4  IVIEW2, 

2  IVIEW3,  IVIEW4 
IVlEUl-0 

IVIEW2-0 
IVIFW3  0 
IVIEW4-0 
30  CONlINUE 

EN»  DO  LOOP 

ERROR_MONO(1)-NhOINTS 
ERROR  MIDI < 1 >«NPOINTS 
QUAH ( 1 ) -NhOINTS 
TO I _RNO < 1 ) -NPOINTS 

PLOT  RESULTS 

IFUPLOT.  EQ.  O)  CO  TO  99 

CAf  L  GRAPH  1  (TOT  RNQ,  ERROR  MONO.  LA8X,  LA8Y*  TITLE#  2#  XMIN.  XMAX,  YMIN,  YMAX  ) 
CAIL  CRAPH2<TGT  RNO#  ERROR  MIDI.  TYPE_LINE,  1  > 

URTTE<6,  1019)  XMIN# XMAX.  YMIN.  YMAX 

I 

CAl  L  0RAPH1  (TOT  RNC.  QUAD.  LAP X «  LABY1.  TITLE#  1#  XMIN,  XMAX,  YMIN.  YMAX) 

UR 1 1 E ( 6,  1019)  XMIN.  XMAX.  YMIN.  YMAX 
99  CONlINUE 

FGRilAT  STATEMEN1S 

1000  FGRHATUX.  'MONO  ERROR'.  14) 

1001  FGRilAT ( IX,  'POSITIVE  MONO  CAL  ERROR  AT  I-',  14) 

1002  FGRilAT  <  '  1  ' ,  AQO ) 

1003  FGRilAT (IX.  '  TOT  RNO  TOT  ANO  MR  ERROR_MONO  ERROR_MIDI 

1  QUAD  IVIEWl  IVIEW2  IVIEU3  IVIEW4') 

1004  FGR.IATC  IX.  4(2X.  F10.  2).2X.  F10.  9,  4<3X,  14,  3X)> 

1009  FGRilAT  (211) 

1006  FGRilAT  <1  OF  10.  3) 

1007  FGRMAT<F10.  2) 

1008  FORMATUX,  'NfcCATIVE  MONOjCAl  ERROR  AT  I-',  14) 

1009  FORMAT (IX, 'IVIEWO  ERROR', 14) 

1010  FGRI1AT<  IX,  'MIDI  ERROR',  14) 

1011  FORMATUX,  2<2X,  F10.  9) ) 

1012  FORMATUX,  'CAL  MONO  ELANG  MR') 

1013  FGRilATUX,  'CAL  MIDI  ELANG_MR ' ) 

1019  FGRilATUX,  'XMIN,  XMAX.  YMIN,  YMAX  '.  4E12.  4) 

1016  FGRilATUX,  'DUM  ERROR  AT  I-',  14) 

1020  FORMATUX.  'MONO  IN_PHASE,  TOT  ANO_MES ',  2<2X,  E12.  9) ) 

1021  FORMATUX,'  C,  TOT  ANG_MES ' ,  2<2X,  E12.  9)  ) 

1023  FGRMAT( IX,  3(2X,  E10.  3) ) 

1070  FORMATUX,  'VOLT  ERRORS',  214) 

1080  FGRMATUX,  'I VIEW  ERRORS  162', 314) 

1081  FGRilAT <  1 X ,  'IVIEW  ERRORS  364',  314) 

1100  FORMAT <A80) 


STOP 

EN» 


tea 

.■  /i 
•*.  *j 


\  ,S 


r.  •  •  »  . 


C 

c 
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non  non  onn  n  n  n  n  n  n 


I  PROGRAM  HULTIDATA  | 

COWARE  CONVENI IQNAL  MONOPULSE  ANO  MIDI 

DIMENSION  SUM  DB(40).D£L  DB(40).ELANQ  DEQ<40). 

1  CAL  MONO (40)* 

2  ELAN0_MR  (  40  )  *  CAL _/1IDI  (  40  )  * 

3  SUM— VOLT <  40  > ,  DEL_ VOLT ( 40 ) 

DIMENSION  QUAD (200) 

DIMENSION  TOT _J?NO ( 200 > »  ERROR _MONO< 200).  ERROR J1IDI  (200) 

CMARACTER*80  TITLE 
CHARACTERS©  LABX 
CH/:RACTER-*40  LABY 
CHARACTERS©  LARY1 
CHARACTERS  TYKE_LIN£ 

COMPLEX  SUMJSIGJDIR*  SUMJBI0J*5ULT»  DEL  SIG  DIR.  DEL  SIQJ1ULT 

COMPLEX  RHQ_CPLX»  A*  B*  SUM_TOTAL.  DEL_TOTAL 

COMPLEX  DIFFUSEJiEAN  DI*  DIFFUSE  SD  DI* DIFFUSE  SD  II. 

1  SPECUL AR_D I .  SPECULAR_I I , 

2  DIFFUSE_M£AN_II 

COMPLEX  DEL_ROTATEt>,  UNIT_IN.  UNITjQUAD 

REA!  IN-PHASE,  MISS_DIST,  FACET_ANO_ANT 

INTERNAL  DATA 

DATA  SUM_SIO_DIR/ ( 1 .  .  0.  ) / 

DATA  DEL  SI0_DIR/(1.  ,0.  >/ 

DATA  UNIT  IN/  ( 1.  .  0.  )/ 

DATA  UNIT  QUAD/  ( 0.  .  1.  )/ 

LARX-'TARCET  RANCE  —  M' 

LARY- 'ELEVATION  ERROR  —  MR' 

LARY1- 'QUADRATURE  COMPONENT  --  DB ' 

TYPEJ-INp* 'DASH' 

INPUT  FROM  ANTENNA  PATTERN  FILE 


READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
READ  ( 
WRITE 
REAn  <! 
WRIIE 
READ  (| 


1009) 
1006) 
1006) 
1006) 
1007) 
1007) 
1007) 
1007) 
1007) 
1007) 
1007) 
1007) 
1111) 
1100) 
1112) 
L 1 13) 


IPLOT.  IBUO 
(SUMJOSd).  1*1.  40) 
<DEL_DB(I).  1*1.  40) 

( ELANO JDEO ( I ) ,  I -1 ,  40 ) 

PHASE_ANQ_DEO 

RHO 

T0T_ALT 
TGT_MAX_RNG 
TOT  MIN  RNC 
RNQ  STEP 
RDR_ALT 
MISSED I ST 

TITLE 

Al  TDIFF 


CONVERSIONS,  ETC 

SQRT2-SQRT  <  2.  ) 
PISSTANd.  ) 
PI_0VhR_2-PI/2. 
DEU_TO_KR*P I  /.  ISO 
SICN—1. 


•V*\- 

'.■■'/•'[••S' 


, 


.. 

l>  .  •  mm_\* 
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non 


DO  9  1-1.40 

IF<Di-LJDB<I>.  LT.  0.  >  SION-1. 

SUN_V0LT<  I  >-10.  ** < SUM J)B <  I ) /20.  > 

DfeL_V0LT< I >-SIGN*10.  *»(DEL_DD< I >/20.  > 

El  AN0_MR  ( I  >-EI  AN0_DE0 <  I  >*DEQ  T0_MR 

9  CONTINUE 

MONOPULSE  CALIBRATION  CURVE 
DO  10  1-1.40 

CAL_MONO< I >-DEL_VOLTd  >/SUM_VOLT< I > 

10  CONTINUE 
C 

C  PREVENT  CALIBRATION  CURVE  FROM  BEING  DOUBLE  VALUED 

C 

DO  19  1-1.39 

IF<CAL_MQNO< I ).  LT.  0.  )  GO  TO  19 
IF<CALJKONO< 1+1 ).  OT.  CAL  MONO < I > )  00  TO  19 
CAL  J10N0d+l  >  * CAL_MONO<  I  >  +.  0001 
WR1TE<6.  1001 )  I 

19  CONTINUE 

DO  16  1-40.  2.-1 

IF(CALJMONO<  I ) .  OT.  0.  >  CO  TO  16 
IF<CALJ<ONO<  1-1).  LT.  CAL_|10N0<  I ) )  00  TO  16 
CAL  MQN0<I-1>*CAL  MONO < I > 000 1 
WRiTE<6. 1008)  I 
16  CONTINUE 

IF(IBUO.  LT.  1)  CO  TO  29 

WRITE <7.  1012) 

DO  22  1-1.40 

WRITE <7,  1011)  C AL_MONO < I ) .  ELANO_MR ( I > 

32  CONTINUE 

29  CONTINUE 
C 

C  MIDI  CALIBRATION  CURVE 

C 

DO  20  1-1.40 

CAL  MIDI  ( I  ) -2.  *-ATAN <  CAL  MONO <  I  > *SQRT2 > 

20  CONTINUE 

IFdBUG.  LT.  1)  CO  TO  24 

WRITE (7,  1013) 

DO  23  1-1.40 

WRlTE<7.  1011 )  CAL  MIDI  d  > .  ELANOJtR  d  > 

23  CONTINUE 

24  CONTINUE 

NP0INT8-0 

C 

C  TABULAR  DATA  HEADERS 

C 

WRITE  <20.  1002)  TITLE 
WRITE  <20.  1003) 

C 

C  BEGIN  ENORMOUS  ITERATIVE  LOOP 

C 

30  CONI INUE 
C 

C  VECTOR  SUMMATION  OF  RETURNS  FROM  FACETS 

C 

C  READ  TARCET  POSITION.  DIRECT  SIGNAL. SPECULAR  RETURN 

C 

READ (30.  ■*>  TOTJ«»  TOT_Y,  TOT_Z,  VEL_TOT_X.  VEL  TOT  Y,  VEL  TOT  Z 


noooo  o  ooo 


.  ■  V-  '  1  I  .  ■  l  .  <1  t  I  .  .  .  »  .1  .  ■■  tf'  .MM 

WRITE <50,  1200) TOT  X,T©T  Y,T©T_Z.  V£l  TCTr_X,  YEL_tOT_Y,  VhL_IO l_Z 
RFAD(30,  *)DIR  DIRjSIQ 
WRlTE<50, 1220)DIR_DIR_SI0 

RF*D<30, *)  SPEC  RNG,  SPEC_ANG,  SPECULAR  DI, SPECULAR _I I 
WRITE <50,  1210) SPEC _RNG,  SPEC_AN0,  SPECULAR_DI,  SPECULAR_II 

1200  FORMAT (IX#  "TOT  X.  TOT_Y,  TOT_Z,  VtL_TGT_X,  VEL  TQT_Y»  VEL_TQT__Z- ' , 

1  6E12.  4) 

1210  FORMAT (IX*  "SPEC  J*NG, SPEC _ANO,  SPECULAR  DI.  SPECULAR_II- ", 

1  6E12.  4) 

1220  FORMAT (IX,  "DIRECT  SIGNAL  VOLT-".  E12.  4> 

TO I  J*NQ (NHOINTS+2 ) -T0T_Y-88. 

TC1_^N0— (ALT_J)IFF)*1000.  /T0T_RN0(NP0INTS+2) 

SPE«:_ANO_ANT— 1000.  *8PEC_AN0-T0T_AN0 

AM)  ANC_DIRECT— 0. 

CAIL  LIN1NT  (ELANO  MR,  SUM_VOLT,  40,  ANT  ANO  DIRECT,  GAIN_SUM» 

1  IVIEWO) 

CAI  L  LIN I NT (ELANO  MR,  DEL  VOLT,  40.  ANT  ANG  DIRECT, GAIN_DEL, 

1  IVIEWO) 

IF<  IVIEWO.  Nfc.  0.  )  WRITE <6.  1009)  IVIEWO 


CALL  LINIMT(E)  ANO  MR,  SUM  VOLT. 40. SPEC  ANO_ANT,  VOLT_MULT_SUM. 

1  I view!) 

SUM  TOT  AL-SPECUL  AR  _D  I  *0  A I  N_SUM  * VOLT_MULT JSUM+SPECUL AR  _I  I 
1  *-VOLT j:ULT_SUM**2+SUM  810  DIR*DIR_DIR_SIO*OAIN  SUM**2 

C,V  L  LIMIMT(E)  ANO  MR,  DEL_VOLT,  40,  SPEC_ANO_ANT,  VOLT_MULT_DEL, 

1  IVIEW1) 

DEL  TOTAL-SPECULAR  D I  *G  A I  N  JBUM  >■  VOLT  MULT  DEL/2.  +SPECULAR  II 

1  *VOLT  hULT_SUM*VOLT_MULT  DEL+DEL_.SIO_DIR*DIRJDIR JBIG 

2  *OAIN  SUM*© A I N JDEL+SPECUL AR  DI*VOLT  MULT_SUM 

3  *0A  I NJDEL  /  2. 

RFAi>  DIFFUSE,  WEIGHT  WITH  ANTfcNNA  OAIN 

READ  <  30.  * )  NFACETS 
WRITE <50. 1230)  NFACETS 
1230  FOPMATUX,  'NFACETS-  ".15) 

DO  35  J-l. NFACETS 

RFAD(30.  *)  RNG  FAC.  EL _ANO_FACET,  AZ_ANO_FACET, 

1  DIFFUSE  ME AN _D I,  DIFFUSE  SD  DI, 

2  DIFfUSE_MEAN_II.DIFFUSE_SD_II 
TWT_FACET_ANG  EL_ANQ_FACET-TCT_ANQ/ 1000. 

FACET  ANO  ANT— 1000.  *ACOS< COS < TOT _pACET  ANO)*COS< AZ_ANO  FACET)) 
CALL  LININTCEI  ANO  MR,  SUM_VCLT,  40,  FACET_ANO_ANT,  VOLT_MULT_SUM, 

1  I VIEW!) 

SUM  TOTAL-SUM  TOTAL+DIFFUSE_SD  D I *CA I N_SUM#VOLT_MULT_SUM 
1  +D IFF USE  SD_I I *VOLT_MULT  SUM**2 

CALL  LININT(E)  ANO  MR,  DEL_VOLT,  40,  FACET_ANO_ANT,  VOLT_MULT_DEL, 

1  IVIEW1) 

DtL  TOTAL-DEL  TOTAL+DIFFUSE  SD  DI*OAIN  SUM#VOLT  MULT_DEL/2. 

1  +D IFF USE  SD  I I *VOLT  MULT  SUM*VOLT_MULT  DEL+ 

2  DIFFUSE  SdJ)I*V0LT_MULT_SUM*GAIN_DEL/2. 

35  CONTINUE 

SUM  TOTAL  AND  CEL  TOTAL  SIGNAL  COMPUTATION  NOW  COMPLETE 
COMPUTE  IN  PHASE  AND  QUADRATURE  COMPONENTS  OF 
D*-L/SUM  RATIO 

SI  CABS (SUM  TOTAL) 

WRITE <50, 1240)  SI  169 


A  A- 


r  1240 


c 

C 

C 


1 


c 

c 

c 


c 

27 


26 


3 
1 

4 
2 


FORMAT* 1*,  'SI-  V6E12.  4) 

Di -CARS (DEL  TOTAL) 

Rl-REAI  (SUM  TOTAL* CON JO ( DEL_TOTAL ) ) 

IN_PHASE-R 1 / (S1++2 > 

AKC-AC0S(R1/(D1*S1)  ) 

QUAD ( NhO I NTS+2 ) *  D1*SIN( ANO) /SI 
QW.H  (  NHQINTS42  )  -QUAD  <  NPOINTS+2 ) «  *2 
IF(QUAn(NH0INTS+2).  LT.  .  0001 )  QUAD (NPOINTS+2)-.  0001 
Q'J/.n(NH0INTS+2)-10.  *AL0Q10 ( GUAM ( NPOINTS+2 )  ) 

MONOPULSE  MEASUHEMFN  T 

CAI  L  L IN INT (  CAL  MONO >  ELANO_MR,  40/  IN_PHASE* 

AN I  TOT  ANOJMES, IVIEW3) 

TCT_ANC  HP  S-AN T_TOT_ANO_MES 

IF(  IBUO.  EQ.  1 )  WHITE (7.  1020)  IN_PHASE,  TOT_ANO_MES 
ERROR  _MONO ( NPO I NTS+2 ) «TQT_AN©  J1ES 

MIDI  HPASUREMENl 

DELJ10TATED-DELJT0TAL/SQRT2 

CAI  L  ROTATE (DEL  ROTATED*  PI JDVERJ2) 

A— SUM__TOT AL/2.  +DEL  ROTATED 
8— SUM_TOT AL/2.  -DEL_ROTATED 
DUM1-REAI  (A-*CONJC(B )  ) 

DUMP-CABS ( A ) *C  ARS ( B ) 

IF(nUr11.  CE.  DUMP)  00  TO  27 
C— ACOS  (  DUM1  /DUM2  ) 

CO  TO  26 

USE  LAST  VALUE  OF  C 
CONIINUE 

WRI1E(6, 1016)  NKJINTS+1 

IF(IBUO.  EQ.  1)  WKITE(7,  1016)  NKJINTS+1 

CAI  l  LININT(CAL  MIDI.  ELANOJMR* 40*  C.  ANT_TOT_ANO_MtS* IVIEW4) 

TCI  ANC_M+  S— AN X  TOT_ANO JMES-ANO _OFF_BORES I GHT 

IF( ( IVIEW3+IVIEW4).  OT.  1 )  WRITE<6,  1081 >  IVIEW3. IVIEW4. NPOINTS+1 

ERROR_MIDI  ( NPOINTS+2 )  — T OT_ANO JflES 

IF (  IBUO.  EQ.  1 )  WRITE(7,  1021 )  C.  TOT_ANO_MES 

NPOINfS— NHOINTS+1 

WRITE (20,  1004)  TOT  RNO ( NPOINTS+1 ) »  TOT_ANO* 

ERROR_MONO( NPOINTS+1 ) . 

ERROR  MIDI (NPOINTS+1 ),  QUAD ( NPOINTS+1 ) . 
IVIEW2. 

I VIEWS, IVIEW4 

IVlEWl— 0 
IVIEW2-0 
IVIFW3  0 
IVIEW4-0 

REAM (30*  ■*)  IFLAO 
WRITE (30,*)  IFLAO 
IF( IFLAO.  EQ.  1)  CO  TO  30 
IF(NPOINIS.  EQ.  1)  CO  TO  99 


C  Eid>  DO  LOOP 

C 

ERROR _MONO( 1 )-NPOINTS 
EPP.UR  MIDKD-NPOINTS 
QU'.n  ( 1 )  -NhOINTS 
TC \  _RNQ ( 1 > -NPO I NTS 

PLOT  RESULTS 


IVlEWl. 


c 

'  c 
c 


IFdPLOT.  EQ.  0)  CO  TO  99 
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1000 

1001 

1002 

1003 

1 

1004 

1005 

1006 
1007 
1006 

1009 

1010 
1011 
1012 
1013 

1015 

1016 
1020 
1021 
1023 
1070 
1080 
1081 
1111 
1100 
1112 

1 

1113 


CALL  GRAPH1  < TOT JW©»  ERROR  MIDI.  LABX,  LABY,  TITLE/  2,  XMIN/  XMAX,  YMIN,  YMAX) 
CALL  0RAPH2(T0T_RN0»  ERROR  MONO.  TYPE  LINE/1) 

WRITE (6* 1015)  XhlNz XMAX.  YMIN,  YMAX 

CALL  GRAPH 1 (TOT  RNO.  QUAD,  LABX,  LABY1,  TITLE,  1,  XMIN,  XMAX,  YMIN. YMAX) 

WRITE (6.  1015)  XMIN,  XMAX.  YMIN.  YMAX 
CONTINUE 

FORMAT  STATEMENTS 
FORMAT (IX.  'MONO  ERROR'. 14) 

FORMAT < IX.  'POSITIVE  MONOjCAL  ERROR  AT  NPOINTS  14) 

FORMAT  (  '  1 ' .  ASO ) 

FORMAT < IX.  '  TOT  RNO  TOT_ANO_MR  ERROR  MONO  ERROR  MIDI 
QUAD  IVIEW1  IVIEW2  IVIEW3  IVIEW4') 

FORMAT ( IX.  4(2X.  F10.  2).2X.  F10.  5.  4<3X,  14.  3X)) 

FORMAT (211) 

FORMAT < lOriO.  3) 

F0RfJAT(F10.  2) 

FORMAT < IX,  'NEGATIVE  MONO_CAL  ERROR  AT  NPOINTS  14) 

FORMATUX,  'IVItWO  ERROR',  14) 

FORMAT < IX.  'MIDI  ERROR'. 14) 

FORMAT (IX,  2(2X,  F10.  5)) 

FORMAT (IX,  'CAL_MONO  ELANG_MR ' ) 

FORMATUX.  'CAL  MIDI  ELANQ_MR  ' ) 

FCRflAT <  IX,  'XMIN.  XMAX.  YMIN,  YMAX  ',  4E12.  4) 

FORMAT UX,  'DUM  ERROR  AT  NPOINTS  «',I4) 

FORMATUX.  'MONO  IN  PHASE,  TOT  AN0J1ES '.  2<2X,  E12.  5)  ) 

FORMAT (  IX.'  C,  TOT  ANGJ1ES',  2<2X.  E12.  5)  ) 

FORMAT <  IX,  3(2X,  E10.  3)  ) 

FORMATUX.  'VOLT  ERRORS',  214) 

FORMATUX.  '  I  VIEW  ERRORS  1&2',3I4) 

FORMATUX.  '  I  VIEW  ERRORS  3*4',  314) 

FORMATUX.  'ENTER  TITLE:  ') 

FORMAT ( ASO) 

FORMATUX.  'ENTER  DIFFERENCE  BETWEEN  TARGET  AND  RADAR  ALTITUDES 
<  fOT  ALT-RDR  ALT)') 

FORMAT <F10.  2) 


•.  •  -»T. ' 


SUBROUTINE  ROTATE  (  XCOMPLEX >  ANGLE  ) 
DIMENSION  XCOMPLEX (2) 

A-XCOM^LEX  < 1 ) 

b*:<cokhlex<2) 

C -COS (ANGLE) 

S"S IN (ANGLE) 

XCOilPLFX  ( 1  )-A*C+B*S 
XCOMPLEX  ( 2 )  —A+S+B-^C 
RETURN 
EN» 


« 


■'IP  A  V.7V.\  '  "  -V"..  - 


I  PP08RAM  CORREL ATI OND I STANCE  VERSION  2  I 

dimension  diffodTcEL  yiSU, HlSSlSi)"  '* 

CHAR AC TER *80  TITLE 
CHARACTER**!)  LADX,  LABY 
CHARACTER**  TYPE_LINE 

DOUBLE  PRECISION  DIST. PATH_LENOTH_DIFF 

REAL  MISS 

LARX  -  'MISS  DISTANCE' 

LARV  -  'DELTA  —  H' 

MR 1 TE ! *< 1000) 

1000  FORMAT! IX.  'ENTER  RDR_ALT.  TOT^ALT.  TGTJ1N0,  DEL  X',4F10.  2) 
RFAD!3.  1001  >rdr_alt.  tot  alt.  tot  rno.  del  X 

1001  FORMAT <4F 10.  2) 

WRITE <6,  1002) 

1002  FORMAT! IX.  'ENTER  TITLE') 

READ <5. 1003)  TITLE 

1003  FCRliAT !  A80  ) 

CALL  DELTA !RDR  ALT.  TOT_ALT. TOT_RNO.  PATH  LENOTH  DIFF) 

C  WRITE!*,  *)PATH_LENOTH_DIFF 

DEI  _Y!1>  -  -13. 

MISS! 1 )  ■  -13. 

DO  10  1-1.30 
MISS! 1*1 )— MIS8! I )*1 
DtL_Y! 1*1 )— MISS! 1+1 > 

RNC-TOT_RNO*DEL  X 
ALT-T8T_ALT+DEL  Y!I+1) 

CALL  DELTA !RDR_^LT.  ALT,  RNO. DIST) 

DIfF ! 1+1 )-DIST-PATH  LENOTH  DIFF 
C  WRITE!*.  *)DIFF!I  +  1) 

10  CONTINUE 

MISS! 1 >-30. 

D I FF !  1 )  —30. 

CALL  CRAPHUHISS.  DIFF.  LABX.  LABY.  TITLE,  1.  XMIN,  XMAX, 

2  YMIN.  YMAX) 


STOP 

Ei!» 


SUBROUTINE  DELTA! HR, HT,  RT.  DEL ) 
DOUBLE  PRECISION  A,  B,  0. Rl. R2, DEL 
A  -  !HT  ♦  HR  >**2 
0  -  1HT  -  HR>**2 
8  -  RT«*2 
Rl  -  DS«RT!A*0) 

R2  -  DSORT ! B*0) 

DtL  -  R1-R2 
RE  I  URN 
EMD 
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PROGRAM  KCORELATE 


Dimension  sumjdbmo).  del_ddmo>,  elano_peo<40). 

CAL_M0N0<40) » 

ELANO  MR<40), CAL.MIDI (40) » 

SUM.VOLT <  40) ,  DEL.VOLT ( 40 > 

DIMENSION  QUAD (200), RELATIVE  ERROR (200) 

DIMENSION  ERROR  J1QN0< 200) .  ERROR_MIDI (200) 

DIMENSION  DEL_Y < 200 ) » DIFF < 200 ).  MISS <200) 

CHARACTER*BO  TITLE 
CHARACTER *40  LABX 
CHARACTER *40  LABY 
CHARACTER*40  LABY1 
CHAR AC TER *6  1 YPEJL I NE 

COW»LEX  SUMJBIOJDIR.  SUM  SIO.MULT,  DEL  SIO  DIR. DEL  SI0J1ULT 
COMPLEX  RHOjCPLX,  SUM.TOTAL,  DEL_TOTAL.  DEL  TOTAL  TWT 
CCWLEX  DELJROTATED.  UNIT.IN,  UNIT  QUAD.  SUM  TOTAL  TOT 
COMPLEX  SUMJBIO.MULT.TOT,  DEL  SIO  MULT  TOT 


'  '  \1 


RFAL  IN.PHASE,  MISS_DIST 

DOUBLE  PRECISION  DIST,  PATH _LFNOTH_DIFF 

WAVfcLENOTH-0.  032 
LA*X-'MIS8  DISTANCE  —  M' 

LARY- 'DELTA  —  M' 

LARYl- 'RELATIVE  ERROR  —  MR' 

IMIERN4L  DATA 

DATA  SUM.SIQ.PIR/ ( 1.  ,  0.  )/ 

DATA  DEL  SI0_DIR/<1.  ,0.  )/ 

DATA  UNIT  IN/<1.  ,0.  )/ 

DATA  UNIT.GUAD/  <0.  ,  1.  )/ 


INHUT 

READ 

READ 

READ 

REAR 

READ 

READ 

REAR 

REAR 

REAR 

REAR 

REAR 

REAR 

REAR 


1005)  IPLOT, IBUO 

1006)  <SUM_DB<I),  1-1,  40) 
1006)  < DEI  JDB (  I ) ,  I ■  1 »  40 ) 

1006)  (ELANO  DEQ<I>, 1-1,  40) 

1007)  PHASE_ANO_DEO 
1007)  RHO 

1007)  TOT  ALT 
1007)  TOT  MAX  RNO 
1007)  TGT_MIN  RNO 
1007)  RNO  STEP 
1007)  RDR  ALT 
1007)  MISS.DIST 
1100)  TITLE 


CONVERSIONS,  ETC 

8QR  T2**SGRT  ( 2.  ) 

PI— 4*ATAN( 1.  ) 

PIjDVfcR  2-PI/2. 
rn.O_T0.MR-P  I/.  180 
felCN— l. 

DO  9  1-1.40 

IF<DEL  DB<I).LT.  0.  )  SION-1. 
SU.1.V0LT <  I  )-10.  ** ( SUM  DB  <  I  >  /20.  > 


V.  V‘v*V 


.  .*  ,*  ,* 


D*-L_VOLT< X )-SICN*10.  **<DELJDD< I )/20.  ) 

El  ANG_MR <  X  )*EL ANOJJEO ( X ) *DEG_TQJ1R 

9  CONTINUE 
C 

C  MQN0PUL8E  CALIBRATION  CURVE 

C 

DO  10  1-1*40 

CALJiQNO < I ) -DEL_VOLT  <  X ) /SUM_VGLT < I ) 

10  CONTINUE 
C 

C  PREVENT  CALIBRATION  CURVE  FROM  BEING  DOUBLE  VALUED 

C 

DO  19  1-1*39 

IF<CAL_r.ONO<I>.  LT.  0.  )  00  TO  19 
IF(CAL_f10N0<I+l>.  QT.  CAL_MONO<I))  GO  TO  19 
CAL  MONO(I+1)=CAL_J10NO<I)+.  0001 
WRITE! 6, 1009)  I 

19  CONTINUE 

DO  16  1—40*  2* -1 

IF(CAL_J1QN0(  I ).  OT.  0.  )  GO  TO  16 
IF<CAL_T.0N0<I-1).  LT.  CAL_NONO(I>>  GO  TO  16 
CA.L  MONO<  1-1  >-CAL_MONO<  I )-.  0001 
UR1TE(6* 1008)  I 
16  CONTINUE 

IF<  I  BUG.  LT.  1)  CO  TO  29 
WRITEC7* 1012) 

DO  22  1-1.40 

WR1TE<7.  1011)  CAL  MONO <  I  ) *  ELANS  .MR  ( I ) 

22  CONTINUE 

29  CONTINUE 

C 

C  MIDI  CALIBRATION  CURVE 

C 

DO  20  1-1*40 

CAL  MIDI ( I ) -2.  *ATAN < CAL_MONO < I ) #SQRT2 ) 

20  CONTINUE 

IF<  IBUG.  LT.  1)  GO  TO  24 

WR1TE(7* 1013) 

DO  23  1-1.40 

UR1TE<7.  1011)  CAL_MIDI < I ) *  ELANO  MR  <  I ) 

23  CONTINUE 

24  CONTINUE 
C 

UR 1 1 E ( 6* 1000) 

C 

1000  FORMAT (IX*  'ENTER  RDR  ALT,  TOT  A!  T,  T0T_RN0,  DEL  X',4F10.  2) 

READ (  9*  1001  )RDRJM_T.  TGT_ALT*  TOT  RNG,  DEL_X 

1001  FGRI1AT<4F10.  2) 

WRITE (6,  1002) 

1002  FORMAT < IX,  'ENTER  TITLE') 

READ ( 9,  1003) TITLE 

1003  FORMAT (A80) 

CAM.  DELTA C RDR _ALT,  T0T_ALT,  TOT  RNQ,  PATH_LENOTH_DIFF ) 

C  WRITE<6,  *)PATH_LENGTH_DIFF 

PHASE  ANGJ?AO-PI+PATHJ_ENOTH  DIFF*2.  *P I /WAVELENGTH 
C 

C  TARCET  ANOLE  ERROR,  MIDI 

C 

TCT_*N0  DIR-1000.  *( TOT_ALT-RDR _ALT  > /TOT_RNO 
TCI _ANO  MULT— 1000.  * ( TQT_ALT+RDR_ALT ) /TOT_RNG 
AN  I  _ANC  MULT-TOT_ANOJ1ULT-TOT  ANOJ)IR 

CAI  L  LIN I NT (ELANO  MR,  SUM  VOLT,  40, ANT_AN0_MULT,  VOLT_MULT_SUM 

174 


o  o  o 


1 


IVi) 

CALL  LININTtEl  ANO.NR,  DEL.VOLT,  40,  ANT.ANO _MULT.  VOLT.MULT .DEL, 

IV2) 

8UM .8  1 O.MULT.TQT- VOLT.MULT_SUH«RHO*UN I T_I N 
DfcL.S  Z  OJlULT.T  OT  *  VOLT  JMULT  JD6  L*RHO*UN  Z  T.X N 
CALL  ROTATE  <SUH_S  I  QJlULT.TQT, PHASE.ANQ  J?AD ) 

CALL  ROTATE  <  DfcL.S 1 ©JlULT.TOT ,  PHA8E_ANG_RAD> 

BOrfES  I  OHT.ANQ*  0. 

CALL  LININTtEl  ANO.NR,  SUM.VOLT,  40,  SORESIGHT _ANO,  BORES I GHT_0 A  X N.SUM, 
1  IV3) 

CALL  LININTtEl  ANQ.MR,  DEL.VOLT,  40,  BORESIOHT _ANQ,  B0RE8IGHT.QAIN.DEL, 
1  IV4) 

SUft.TOTAL.TOT-SUM.S  I Q.D I R*BORES I  CHT.Q A I  N.SUM+ 

1  SUM.SIO  .fUJLT.TOT 

DfcL.TOT AL.TQT-DEL.S  I Q.D I R *B OR ES I  GHT_0  A I  NJDEL+ 

1  DEL.S I Q.MULT.T GT 

SI  CABStSUM.TOTAL  TOT) 

01  CABS  <  D£L_TOTAL_TOT ) 

R1  RFAL  <  SUM  TOTAL  TOT*CONJO ( DEL.TOTAL.TQT ) > 

C  Rl/Si*>2 

CAU  LININTtCAL  MIDI,  ELANO.MR,  40,  C,  ANT  TQT.ANGJ1ES, 

1  IV3) 

TCT_EHROR-AN  r  .TCT  ANO.MES 
WRIlEt023,  1301)  TOT  ERROR,  C,  PHASE_ANG_RAD 

1301  FORMAT <  IX,  'TOT  ERROR,  C,  PHASE.ANO  RAD-',  7F10.  4) 

WRITE <023,  1302) 

1302  FORMAT ( IX,  'DEL  Y  ANT  ANOJDIR  GAIN.DEL 

1  C  AN  I  PROJ_ANC  h>S  PHASE  ANO  RAD') 

DEI  .Yti)--13. 

DO  30  1-1.30 
DhL.Yt 1+1 >-DEL_Y< I )+l. 

RNC-TOT  RNQ+DEL.X 
ALT-TOT  ALT+DEL.Ytl+l) 

CALL  DELTA <RDR_ALT,  ALT,  RNC.  DIST) 

PRGJECTILE  ANGLE  ERROR.  MIDI 

PHASE  ANO  RAD=PI+DIST*2.  *P I /WAVELENGTH 
ANO JDTFJORESICHT— 1000.  *<DEL_Y<I+1  )-DEL_X* 

1  T0T_AN0_DIR/1000.  )/RNQ 

PROU  ANO  DIR-1000.  *-tALT-RDRj*LT)/RNO 
PPOJ  ANO  MULT— 1 000.  * <  ALT+RDR.ALT )  /RNO 
ANT^ANOJD I R»ANO_OFF_BORES I GHT 

ANT.ANG.MUL  T-PR  0 J.ANQ.MUL  T -P R  0 J.ANO.D  I  R+ANC.OFF.BORESIGHT 

CAU  LININTtELANO  MR,  SUM  VOLT,  40.  ANT  ANO.MULT ,  VOLT  MULT  SUM , 

IVI) 

CAU  LININTtELANO  MR,  DEL  VOLT.  40,  ANT  ANO J1ULT,  VOLT  MULT  DEL, 

IV2) 

8UM.8 1 0_MULT  -VOLT.MULT.SUM-aRHO^UN  I  T.I  N 
DEI  _S I C_MULT-V0LT_MULT_DEL*RH0-*UN I T_I  N 
CAI  L  ROTATEtSUM.SIO.MULT,PHASEJWQ.RAD) 

CAI  L  ROTATE  <  DEL.S  I  C.MVLT,  PHASE.ANQ.R AD ) 

CAI  L  LININTtELANO.MR,  SUM.VOLT,  40,  ANT  ANOJDIR, 

QA I  N.SUM.  IV3) 

CAI  L  LININTtELANO  MR,  DEL.VOLT,  40,  ANT.ANQ.DIR, 

QAINJDEL, IV4) 

SUM.TOTAL  -SUM  S 1 0_D  I  R*OA  I  N.SUH+SUM.S  1 0_MULT 
DEI  .TOTAL -DEL  S I C_D I R*OA I N.DEL+DEL.S 1 0_MULT 
8 1 -CARS  t SUM.TOT  AL ) 

D1-CAB8 1  DEL  .TOTAL ) 
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R 1  RF  AL  <  8UM_TQT AL»C0N0Q < DEL_TOT AL ) ) 

C‘R1/81*»2 

CAI  L  UNIMTtCAL  J1IDI.  ELANO_MR,  40.  C.  ANT_PRCW _ANQ  J1ES, 

1  IV6) 

PROJ_EKROR-ANTJ>ROUJWQJ1ES-ANQJ3FFJBORE8IOHT 
REI  ATIVEJERROR ( X+l >  »PROU_ERROR-TQT_ERROR 

WRI  IE<025,  1300)  DEL_Y(I+1 ) ,  ANTJANO.DIR,  GAIN  DEL.  C,  ANT^PROJ _jANO  MES 
1  .  PHASE_jANO_RAD 

1300  FORMAT (10F10.  4) 

DIFf"<  1+1 )-DIST~PATHJ.ENQTH_DJFF 
C  URiTE(6,*)DIFF(I+l) 

30  CONTINUE 

DEI  _Y(  1  )»30. 

DIFF(1>«30. 

REi  ATIVEJERROR  ( 1 ) -30. 

CALL  CRAPH1  (DfcLJV,  DIFF,  LABX.  I  ABY,  TITLE.  1,  XMIN,  XMAX. 

1  YM1N.YMAX) 

CALL  GRAPH  1  (DfcLJY,  RELATIVE^. ERROR*  LABX,  LABY1,  TITLE,  1,  XMIN, 

1  XMAX,  YHIN,  YMAX) 

FORMAT  STATEMENTS 

1005  FORMAT <21 1) 

1006  FORMAT < 10T10.  3) 

1007  F0RMAT(F10. 2) 

100B  FORMAT! IX,  'NEGATIVE  MONOjCAL  ERROR  AT  I-M4) 

1009  FORMAT < IX,  'POSITIVE  MONO  CAL  ERROR  AT  I-',  14) 

1011  FORMAT < IX,  2(2X,  F10.  5) ) 

1012  FORMAT (IX,  'CAL_MO.NO  ELANQ  MR') 

1013  FORMATOX,  'CAL_MIDI  El  AMO  MR ' ) 

1100  FORMAT (ASO) 

STCr 

END 

SUHROU!  IMS  DELTA  (HR,  Hr,  RT,  DEL) 

DOUBLE  PRECISION  A, B, G, Rl, R2, DEL 
A  «  (HI  +  HRH»2 
B  -  (HI  -  HR ) -**2 
G  -  RT**2 
Rl  »  DSGRT ( A+C ) 

R2  -  DSGRT CB+C) 

DO  -  R1-R2 
RE I URN 
EM» 


C 

SUDROUT IKE  ROTATE ( XCOMPL EX .  ANCLE ) 
DIMENSION  XCOMPLEX ( 2 ) 

A^XCOMPLEX(l) 

B» XCOMPLEX (2) 

C  COS (ANGLE) 

S-S IN (ANCLE) 

XCOMPLEX <1)«  A«C+B*8 
XCOMPLEX  ( 2 )  —A*  S+B*C 
RE  I URN  175 

— ^  _  END 
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I  PROGRAM  1NTERFER_T0_8I0  j 

COi<PUT E  CLUTTER  TO  TARGET  AND  MULTIPATH  TO  TARGET  RATIOS 

DIMENSION  SUM_DB< 16) <  DEL  DB( 16) >  EL ANO_DEQ < 1 6 > . 

1  ELAND  MR ( 16) 

DIMENSION  TQTJTN0<200).  MULT_TO_TGT_SUM < 200 ) . 

1  CLUT_TO_TQT_SUM <  200 ) .  MULT_TO_TGT_DEL <  200  > . 

2  CLUT_TO_TQT_DEL  <  200  > 

CHARACTFR*80  TITLE 
CHARACTER*40  LABX 
CHAR ACTE  R*40  LABY 
CHARACTER*40  LABY1 
CHARACTERS  IYPEJLINE 

REAL  MULT_TO_TQT_SUM. MI83_DIST,  MULT_T0_T0T_DEL 

INlERNAL  DATA 

LAKX-'TAROET  RANGE  —  M' 

LABY- 'CLUTTER  TO  TARGET  RATIO  -DB ' 

LADYl-'MULTIPAIH  TO  TARGET  RATIO — DB ' 

TYPE  J..INE—  'DASH ' 


INr'l/T  FROM  FILE 


RFAD  (10.  1009) 
RFAB  <10.1006) 
RFAD  <10.1006) 
RFAD  <10.1006) 
RFAD  <10.1007) 
RFAD  <10.1007) 
RFAD  <10.1007) 
RFAD  <10.1007) 
RFAD  <10.1007) 
RFAD  <10.1007) 
RFAD  <10.1007) 
RFAD  <10.1007) 
RFAD  <10.1100) 
WRITES.  lOOO) 
RFAD  <  S. 1001 ) 


1PL0T.  I BUG 
<8UH_DB< I >.  1-1.  16) 
<DEL_DB( 1 ).  1-1.  16) 
<bLANGJDEO< I ). 1-1.  16) 
PHASE_ANG_DEG 
RHO 

TOT_ALT 

TOT_MAX_RNO 

T0TJ1INJTN0 

RN0_8TEP 

RDR_ALT 

MISSJJIST 

TITLE 

TITLE 


INPUT  FROM  TERMINAL 
URlTE<6. 1200) 

RFADO.  *>SI8MA_ZER0.  8IGMA.TGT.  BEAM  WIDTH.  RNGJIES 
FORMAT (IX,  'ENTER  SICMAJtEPO— DD,  8I0MA_T0T~ DBSM, 
BFAM_WI  DTH-  -DEO.  RNO_RES — f1 ' ) 

WRITE) 10.  * > 8 1 CMA_ZERO.  810MA_TGT.  BEAM  WIDTH, RNGJIES 

COWERS  IONS.  ETC 

SA|<T2-S0RT<2.  ) 

P I -4* AT AN ( 1 .  > 

PI  0VFR_2-PI/2. 

D*-3  TO  MR -PI/.  180 
DO  9  1-1.  16 

El  /NO  NR  < I )-£(  ANO  DEO < I ) *DEG_TO_MR 
CONTINUE 

DF.\?  i_U  I DT  H_R  AD-DE  AM_W  I DT  H*DCO_TO_MR  / 1 000 . 

NPOINTS  0 


•  •  •  *»  . 
f  V  v 

*  \*  ol 
•• 
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1 

1 


1 

1 

1 

1 

1 

1 

1 

30 


1 

1 

99 

1000 


TCTJINO! 1 >-T0T J1IN  RN0-RN0_STfcP 

MSTfcPS-INI  ((TCTJIAX  RNO-TOTJ1IMJINQ>/RNO  8TEP1+1 

IF(  IVIEUO.  NF.  0.  >  WRITE <6.  1009)  IVXEW0 

BECIM  CALCULATION  00  LOOP 

00  30  I-l.NSTEPS 

TCI  JIN0(  I+l  )-T0TJ»N0<  I  >+RW0_STfcP 

ANS  OFl-  BORES X OHT- 1000.  PMISS  J)IST/TQTJINQ(  I+l  > 

AN I  _ANC  DIRECT«>ANO_OFF JJORESICHT 

CAI  L  DQUAO I N T < O  ANC_HR .  SUM  DB.  16.  ANT_ANO  DIRECT.  ©AIN_SUM, 

IVIEUO) 

CAI  1  DQUAOINT (El  ANOJfR.  DEL_DS.  16.  ANT JMWJJIRECT,  OAXN_DEL. 

IVIEWOO) 

T0TJkNQJ)IR-1000.  * < TOT_ALT-RDR JILT+M I  SB  J)  I  ST  )  /  TOT  J!  NO  <  I + 1 ) 
TQT_AN0J1ULT— lOOO.  *  ( TOT_ALT+RDR  JILT+M  X  88  JD 1 8T  )  /TOT_R  NO  <  I+l  > 

ANT  _AN3  JflJLT-TOT JIN0J1ULT+AN0_0F F  J90REB I CHT-TOT_ANO  JD  X  R 
Cl  UT _ANC— 1000.  +RDR_^LT /TOT  JINS  (X+l) 

ANT_AN3_CLUT-CLUT_AN0-T0T _ANO  DIR+ANO_OFF  DORESIQHT 
CALL  DQUAOINT  (ELAKO  |iR.  SUMJ)D.  16.  ANT  JINO  JKJLT,  ©AINJHJLT. 

IVIEU1 ) 

CALL  DGUADINTIELANOJIR.  8UM_DD.  16.  ANT _ANO  CLUT.  OAIN_CLUT, 

IVIEU2) 

8JCMA_CLUT-10.  +AL0310I TOT_RNO< 1+1 >*BEAM  WIDTH JIAD+RNO  JIEB) 
+81CNA_ZER0 

CLUT_TO  TOT_BUH( X+l )«2+(0AIN  CLUT-8AXN_SUM>-8ICMA_T0T+SI0MAjCLUT 
MliLT_TO  TQT_BUH(I+1)»QAINJ1ULT-QAIN  8UM+20.  •ALOOIO(RHO) 

CALL  DQUAOINT (ELANS  JtR.  OELJDD.  16.  ANT_ANO  J1ULT .  SAIN_MULT_DEL, 
IVIEWO) 

CALL  DQUAOINT  (ELANQJIR.  DEL_DB.  16. ANT_AN0J1ULT.  SAIN_CLUT  DEL. 
IVIEU4) 

Cl  UT  TO  TOT  D+L(I+1)»0AIN_CLUT„SUH+0AIN  CLUT  DEL-OAIN  DEL 
-SIQMA_TOT+SIOMA  CLUT 

Ml.-LT  TO  TOT  D+L(  1+1  >»OAINJlULT  DEL-OAIN  DEL+20.  *AL0010(RH0> 

WRITE! ll.+JTQT  RN0( I+l ).  CLUT  TO_TQT  SUnTl+1 ).  OAIN_SUM. 

OAJNjCLUT.  8I0MAJCLUT 

NP0INT8-NP0INT8+1 

CONTINUE 

END  DO  LOOP 

CLUT_TO_TOT JBUM ( 1 )-NPOINTS 
MULT_TO  T0T_8UM(1>-NP0INT8 
CLUT_T0_TQT_PEL ( 1 J-NPOINTS 
MULT_TO_TOTJ>EL(  1  )-NPOINTS 
TOT  J(NS( 1 )-NP01NT8 

PLOT  RE8ULT8 

C.M  L  GRAPH 1 (TOT JINS.  CLUT_TO_TOT_SUH.  LABX.  LABY.  TITLE.  2.  XMIN,  XMAX. 
YMIN.  YMAX ) 

CALL  CRAPH2( TOT  JINS.  CLUT_TO  TOT  DEL. TYPE  LINE.  1) 

CALL  CRAPH1 (TOT JINS. MULT_TO_TOT_SUM.  LABX.  LABYl.  TITLE,  2.  XMIN.  XMAX. 
YMIN. YMAX) 

CALL  CRAPH2(T0TJIN0.MULT_T0_T0T_DEL.  TYPE  LINE. 1) 

WRITE(6.  1019)  XMIN.  XMAX,  YMIN,  YMAX 

CONTINUE 

FORMAT  STATEMENTS 
FGRilAT (IX,  'ENTER  TITLE:  ') 
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Figure  D-2.  MIDI  tracking  error:  York  scoring;  X-band;  40dB  MTI  improve 
ment  factor;  target  height  *  50  m. 
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Figure  D-3.  MIDI  tracking  error:  York  scoring;  X-band;  50dB  MTI 
improvement  factor;  target  height  =  50  m. 
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RMS  ERROR  (MILS) 

AM  IM  140 _ t-tO  IM  200  2  20  2  40  2.  IK)  K>  .1 i.00 
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*>0°  120  00  ’Woo  200.00  2'^o.oo  2'ao.oo  3io!oo  1mm  4>o.oo  440.00 
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Figure  D-4.  MIDI  tracking  error:  York  scoring;  X-band;  40dB  MTI 
improvement  factor;  target  height  =  15  m. 
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Figure  0-6. 


MIDI  tracking  error:  Missile  scoring;  Ka-band;  40d8  MTI 
iMprovsMant  factor;  target  height  «  150  m. 
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